AD-A047  261 


UNCLASSIFIED 


MICHIGAN  UNI V ANN  ARBOR  DEPT  OF  NAVAL  ARCHITECTURE  — ETC  F/G  20/4 
OPTIMAL  STOCHASTIC  PATH  CONTROL  OF  SURFACE  SHIPS  IN  SHALLOW  WAT— ETC (U) 
AUG  77  MG  PARSONS#  H T CUONG  N00014-76-C-0751 


DEPARTMENT  OF  NAVAL  ARCHITECTURE  AND  MARINE  ENGINEERING 


THE  UNIVERSITY  OF  ^MICHIGAN 
ANN  ARBOR,  MICHIGAN  48109 


Approved  for  public  release;  distribution  unlimited. 


PREPARED  FOR  THE 

OFFICE  OF  NAVAL  RESEARCH*800  N.  QUINCY  ST . • ARLINGTON *VA*  22217 

//cj 


I 


Unclassified 


SECURITY  CLASSIFICATION  OF  this  PACE  fWi rn  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

before  completing  form 

1.  REPORT  NUMBER 

v/ 

ONR-CR215-249-2F 

2.  GOVT  ACCESSION  NO. 

S.  RECIPIENT'S  CAT  ALOG  NUMBER 

«.  TITLE  (and  Subtitle) 

Optimal  Stochastic  Path  Control  of 
in  Shallow  Water 

Surface  Ships 

S.  TYPE  OF  REPORT  & PERIOO  COVERED 

Final  report 

1 May’ 76  - 30  June  77 

! 

6.  PERFORMING  ORG.  REPORT  NUMBER 

No.  189 

7.  AUTHORS 

Michael  G.  Parsons 

Hua  Tu  Cuong 

0.  CONTRACT  OR  GRANT  NUMBERfa) 

N00014-76-C-0751 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  Naval  Architecture  and  Marine  Engr.1^ 
North  Campus  - The  University  of  Michigan 

Ann  Arbor,  MI.  48109 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  0 WORK  UNIT  NUMBERS 

>1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research,  Code  211 

12.  REPORT  DATE 

15  August  1977 

800  North  Quincy  Street 

Arlington,  VA.  22217 

13.  NUMBER  OF  PAGES 

viii  + 75 

<4-  MONITORING  AGENCY  NAME  ft  AODRESSfJf  di  Iterant  front  Controlling  Ollica) 

IS.  SECURITY  CLASS,  (ol  this  report) 

Unclassified 

I5«.  DECLASSIFICATION/ DOWNGRADING 
; SCHEDULE 

t9.  KEY  WOROS  (Continue  on  reverse  side  It  necessery  end  Identity  by  block  number) 

stochastic  control  systems  Kalman-Bucy  filters 

optimal  control 

surface  ship  path  control 

shallow  water  maneuvering 

modern  control  theory  


TRACT  (Continue  on  reverse  side  It  necessery  end  Identity  by  block  number ) 

The  control  of  a surface  ship  along  a prescribed  straight-line  path  is  formu- 
lated as  a stationary,  linear,  state-variable  control  problem.  The  open-loop 
characteristics  of  this  problem  are  studied  using  data  for  a Mariner  type  ship  at 
two  speeds  and  varying  water  depth-to-draft  ratio  and  for  the  tanker  Tokyo  Mar u at 
one  speed  and  varying  water  depth-to-draft  ratio  using  data  obtained  by  Fujino. 
Optimal  stochastic  control  systems  usi-ng  a Kalman-Bucy  filtef  and  a state-feedback 
controller  are  designed  for  both  vessels  at  various  conditions.  These  designs  are 
developed  to  control  the  ship  when  subject  to  random,  zero-mean  disturbances.  The^ 


DD 


FORM 
JAN  73 


EOITION  OF  1 NOV  63  IS  OBSOLETE 


Unclassified  * 


•SECURITY  CLASSIFICATION  OF  THIS  PACE  Dote  Enter* 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  RAOEf Wh,n  Dal*  Enl*r*d) 

“ design  disturbances  are  the  yaw  moment  and  sway  force  due  to  a passing  ship  which 
are  modeled  by  first-order  shaping  filters  in  the  design  derivation.  ..The  selection 
of— a set  of  measurements  adequate  to  provide  effective  path  control — •fq  studied. 

-)System  performance's  studied  by  the  evaluation  of  the  Root  Mean  Square  (RMS) 
response  of  the  controlled  ship  to  the  modeled  design  disturbances  and  by  digital 
computer  simulation  of  the  response  of  the  controlled  ship  to  initial  condition 
errors  and  the  specific  disturbances  due  to  a passing  ship.  The  effects  of  vessel 
speed  and  water  depth  on  the  design  and  performance  of  these  controllers  are 
studied  in  detail.  These  results  yield  guidance  for  the  selection  of  design 
conditions  for  the  design  of  constant -gain  controllers  and  provide  an  assessment 
of  the  need  for  adaptive  controllers  which  can  adjust  the  gains  to  remain  optimal 
as  the  ship  characteristics  change  with  vessel  speed  and  water  depth.  A sensitivity 
study  is  presented  to  show  quantitatively  which  of  the  system  coefficients  which 
change  with  water  depth  can  be  assumed  constant  and  which  of  the  coefficients  must 
be  considered  variables  in  an  adaptive  path  control  system  design.  \ 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PACEfWhin  Dal*  Erw.r.JI 


II 


No.  189 

Report  ONR-CR215-249-2F 
15  August  1977 


I 


Optimal  Stochastic  Path  Control 
of  Surface  Ships  in  Shallow  Water 


Michael  G.  Parsons 
Hua  Tu  Cuong 


I 


This  research  was  carried  out  under 
ONR  Contract  N00014-76-C-0751 


Department  of  Naval  Architecture 
and  Marine  Engineering 
College  of  Engineering 
The  University  of  Michigan 
Ann  Arbor,  Michigan  48109 


J 


Abstract 


The  control  of  a surface  ship  along  a prescribed  straight-line  path  is 
formulated  as  a stationary  linear,  state-variable  control  problem.  The  open- 
loop  characteristics  of  this  problem  are  studied  using  data  for  a Mariner 
type  ship  at  two  speeds  and  varying  water  depth-to-draft  ratio  and  for  the 
tanker  Tokyo  Maru  at  one  speed  and  varying  water  depth-to-draft  ratio  using 
data  obtained  by  Fujino.  Optimal  stochastic  control  systems  using  a Kalman- 
Bucy  filter  and  a state-feedback  controller  are  designed  for  both  vessels 
at  various  conditions.  These  designs  are  developed  to  control  the  ship  when 
subject  to  random,  zero-mean  disturbances.  The  design  disturbances  are  the 

yaw  moment  and  sway  force  due  to  a passing  ship  which  are  modeled  by  first- 

J 

order  shaping  filters  in  the  design  derivation.  The  selection  of  a set  of 
measurements  adequate  to  provide  effective  path  control  is  studied.  System 
performance  is  studied  by  the  evaluation  of  the  Root  Mean  Square  (RMS)  response 
of  the  controlled  ship  to  the  modeled  design  disturbances  and  by  digital 
computer  simulation  of  the  response  of  the  controlled  ship  to  initial  condition 
errors  and  the  specific  disturbances  due  to  a passing  ship.  The  effects  of 
vessel  speed  and  water  depth  on  the  design  and  performance  of  these  controllers 
are  studied  in  detail.  These  results  yield  guidance  for  the  selection  of 
design  conditions  for  the  design  of  constant -gain  controllers  and  provide  an 
assessment  of  the  need  for  adaptive  controllers  which  can  adjust  the  gains 
to  remain  optimal  as  the  ship  characteristics  change  with  vessel  speed  and 
water  depth.  A sensitivity  study  is  presented  to  show  quantitatively  which 
of  the  system  coefficients  which  change  with  water  depth  can  be  assumed 
constant  and  which  of  the  coefficients  must  be  considered  variables  in  an 
adaptive  path  control  system  design. 
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1.  Introduction 

The  problem  of  controlling  surface  ships  along  prescribed  paths  in 
maneuvering  situations  is  becoming  increasingly  important  from  operational, 
safety,  and  environmental  viewpoints.  In  these  situations,  the  ships  are 
subjected  to  short-term  disturbances  due  to  passing  ships,  current  and  wind 
changes,  waves,  and  bank  and  bottom  changes.  To  complicate  the  problem,  the 
dynamic  characteristics  of  surface  ships  change  significantly  with  depth- 
under-keel  and  vessel  speed.  These  maneuvering  situations  can  place  severe 
demands  on  pilots,  conning  officers,  and  helmsmen.  Thus,  some  form  of 
automated  path  control  may  become  highly  desireable  or  necessary  in  the 
future.  An  adaptive  system  which  can  automatically  account  for  the  changes  in 
ship  characteristics  with  changing  depth  and  speed  may  be  needed. 

With  the  increasing  size  of  ships,  more  meaneuvering  is  now  performed 

in  waters  which  must  be  considered  shallow;  i.e.,  when  the  ratio  of  water 

depth  H to  ship  draft  T becomes  less  than  3 or  so  and  significant  changes 
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in  maneuvering  characteristics  occur.  Fujino  ' has  shown  that  as  ships 
enter  shallow  water  they  usually  become  less  course  stable  as  H/T  approaches 
about  2 and  then  become  more  course  stable  in  even  shallower  water.  Some 
ships,  such  as  the  150,000  deadweight  ton  tanker  Tokyo  Maru  tested  by  Fujino  , 
are  course  stable  for  large  and  small  values  of  H/T  but  are  unstable 
for  some  intermediate  range  of  depth-to-draft  ratio.  The  water  depth  would 
then  be  expected  to  have  a significant  effect  on  the  design  and  performance 
of  surface  ship  path  controllers. 

Most  surface  ship  autopilots  or  controllers  in  use  today  are  designed 
to  provide  only  heading  control  in  deep  water.  These  are  usually  analog 
controllers  which  utilize  various  dead  bands  and  operator  adjustments  to 
improve  performance  in  waves,  adapt  for  speed  changes,  etc.  In  the  maneuvering 
situations  of  interest  here  there  is  a need  to  control  along  a prescribed  path 
rather  than  just  maintain  a prescribed  heading.  Similar  needs  exist  in 
underway  replenishment  and  mine  sweeping.  In  path  control,  the  objective 
is  to  control  the  deviation  of  the  ship  from  the  prescribed  path  as  well  as 
the  heading.  A multiple-variable  control  problem  such  as  this  is  best  treated 
using  modern,  state -variable  control  methods.  These  methods  are  directly 
compatible  with  the  use  of  digital  computer  controllers.  Further,  since 
most  of  the  disturbances  which  the  ship  experiences  are  random  and/or  unpre- 
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dictable  in  advance  it  can  be  helpful  to  treat  the  ship  motion  and  disturbances 
as  stochastic  processes  when  developing  the  control  systems. 

The  theory  for  the  optimal  control  of  linear  stochastic  systems  was 

3 4 5 

well  developed  in  the  1960's.  Presentations  by  Bryson  ' and  Gelb  are 
particularly  clear  and  useful.  Timman^  introduced  the  use  of  optimal  stochastic 
control  systems  for  ship  control  in  this  country  in  1972.  Some  investiga- 
tions concerning  the  use  of  stochastic  control  systems  for  surface  ship 
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path  control  have  recently  appeared  outside  the  United  States.  Zuidweg  ' 
developed  an  optimal,  stochastic  path  controller  for  a ship  subjected  to 

g 

random  disturbances  due  to  waves,  wind,  and  current  in  deep  water.  Millers 
developed  an  optimal,  stochastic  path  controller  for  the  maneuvering  of  a 
large  tanker  through  a narrow  passage  involving  a slight  turn  in  the  approach 
to  Brofjorden,  Sweden.  We  feel  the  use  of  optimal,  stochastic  control 
systems  such  as  these  offers  considerable  promise  for  improved  surface  ship 
control  and  safety  in  the  future. 

There  have  been  primarily  three  objectives  in  the  research  which  we  have 
recently  completed. 


• the  formulation  and  general  characterization  of  the  surface 
ship  path  control  problem  in  shallow  water  using  the  state 
variable  approach; 

• the  study  of  the  feasibility  and  effectiveness  of  using  optimal 
stochastic  controllers  for  surface  ship  path  control  in  shallow 
water  including  an  evaluation  of  the  effects  of  speed  and  depth 
changes  on  controller  design  and  performance  and  an  assessment 
of  the  need  for  adaptive  systems;  and 

• the  development  and  careful  documentation  of  computer-aided 
tools  for  the  design  and  evaluation  of  optimal  stochastic 
controllers. 

This  paper  presents  the  results  of  the  first  two  portions  of  this  work. 

Section  2 presents  the  derivation  of  the  linear  state-variable  equations  of 
motions  for  the  movement  of  a surface  ship  along  a straight-line  path. 

The  effect  of  water  depth  and  vessel  speed  on  the  open-loop  characteristics  of 
the  Mariner  type  ship  and  Tokyo  Maru  studied  Fujino  is  also  presented.  In 
Section  3 we  formulate  the  optimal  stochastic  control  problem  and  present 
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the  results  for  the  design  of  a controller  for  the  Tokyo  Mar u at  one  of  its 
least  course  stable  depth-to-draft  ratios.  Section  4 presents  an  investigation 
of  the  effect  of  ship  speed  and  water  depth  on  the  design  and  performance  of 
the  optimal  stochastic  path  controllers.  In  Section  5 we  briefly  discuss  the 
digital  computer  simulation  of  stochastic  controllers  and  present  selected 
simulation  results  for  the  controllers  we  have  designed  for  the  Tokyo  Maru. 
Finally,  in  Section  6 we  present  an  evaluation  of  the  need  to  adapt  path  control 
systems  for  water  depth  dependent  changes  in  the  various  terms  in  the  equations 
of  motion. 

The  final  objective  of  our  research  has  been  the  development  of  computer- 
aided  design  tools.  Much  of  our  work  utilized  the  OPTSYS  program'*' ^ which  was 
originally  developed  under  A.E.  Bryson,  Jr.  at  Stanford  University  by  W.E. 

Hall,  Jr.  and  extended  by  other  Stanford  students.  We  gratefully  acknowledge 
Prof.  Bryson  for  providing  the  version  of  this  program  which  we  have  modified 
and  adapted  for  our  use  under  the  Michigan  Terminal  System  (MTS) . This 
program  provides  very  fast  and  efficient  solutions  to  the  steady-state, 
linear  optimal  control  and  filtering  problems  including  root  mean  square 
(RMS)  response  evaluation  and  controllability  and  observability  evaluation. 

The  principal  new  tool  developed  in  support  of  the  work  reported  here  is  the 
SHIPSIM/OPTSIM  program  for  the  simulation  of  the  response  of  stationary  linear 
optimal  stochastic  control  systems  to  initial  condition  errors  and  specific 
process  disturbances  while  subject  to  measurement  noise.  Our  separate 
report  on  SHIPSIM/OPTSIM'*''*'  includes  user's  documentation,  programmer's 
documentation  and  listings  for  this  program  plus  user's  documentation  for 
our  version  of  the  OPTSYS  program.  This  reference  is  a useful  supplement 
to  the  material  presented  here. 


2 . Equations  of  Motion 

The  development  of  the  linearized,  state-variable  equations  of  motion 

for  a surface  ship  moving  in  the  horizontal  plane  presented  here  begins  with 

2 

the  formulation  presented  by  Fujino  except  that  a canal  wall  is  omitted. 

The  coordinate  system  for  the  problem  is  shown  in  Fig.  1.  The  O-Ch  system 
is  fixed  in  space  with  the  desired  straight  ship  path  along  the  £-axis. 

The  CG-xy  system  is  fixed  at  the  center  of  gravity  of  the  ship.  The  direction 
of  the  arrows  define  the  positive  sense  of  the  drift  angle  6,  the  heading 
angle  i|i,  and  the  rudder  angle  6.  We  consider  here  only  motion  in  the 
horizontal  plane  and  assume  that  roll  is  negligible.  Thus  only  three 
coordinates  x,  y,  and  <p  need  to  be  considered. 


U A 6 / x , X , u 


y,Y,v 


u = dx/dt 
v = dy/dt 

U = (u2+v2)1/2 

r = diji/dt 


Figure  1 . Coordinate  System  for  Path  Control 


The  exact  equations  for  the  horizontal  plane  movement  of  a ship  are 

integro-dif ferential  equations  in  which  convolution  integrals  represent  the 

i o 1 -j 

~ 1- . — C C A.  — C 1 J j 1 4 / TV  1 4-  1 . . 


history  effect  of  the  fluid  to  previous  ship  motion.' 


Alternatively, 


the  equations  can  be  written  as  differential  equations  in  which  the 
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coefficients  are  frequency  dependent.  Recently,  Fujino  has  shown  that 
for  maneuvers  of  interest  here  the  frequency  dependence  is  negligible  and 
constant-coefficient  differential  equations  can  be  utilized.  For  small 
deviations  as  will  be  true  with  effective  path  control,  the  differential 
equations  can  be  linearized  about  the  ni  minal  path  by  taking  small  perturbations 
from  the  equilibrium  state  and  retaining  only  the  linear  terms.  Applying 
the  usual  symmetry  arguments,  these  linearized  equations  can  be  reduced  to  the 
following^ : 
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(m+mx)^-  = Xuu, 

dv  • i 

(m+my)|^  = Yvv  + (-mU+Yr) r + Yf-r+Y66+Yj6 

dr  * 

dzz+Jzz)^  = Nvv  + Nrr+N^v+N66+N^6  , 


(1) 

(2) 

(3) 


where  the  derivative  Xu  is  the  partial  derivative  of  the  surge  force  with 
respect  to  a longitudinal  velocity  perturbation  u with  all  other  quantities 
at  the  equilibrium  straight-line,  constant  speed  U condition,  etc.  In  the 
linearized,  small  perturbation  equations,  the  longitudinal  or  surge  equation 
decouples  and  thus  need  not  be  considered  further  here.  Further,  the  terms 
involving  6 are  small  and  can  be  neglected. 

For  the  path  control  problem,  we  are  also  interested  in  the  lateral  offset 
or  deviation  from  the  desired  path,  n,  which  is  given  by  the  additional 
differential  equation, 


~ = U sin  (i|>-8)  , 
at 


(4) 


for  which  the  following  approximation  is  valid  for  small  deviations, 

= U (*-B) . (5) 

It  is  common  to  utilize  the  drift  angle  instead  of  the  lateral  velocity 
in  ship  maneuvering  work.  These  are  related  by, 

v = -U  sin  8,  (6) 

for  which  the  following  approximation  is  valid  for  small  deviations, 

v = -UB  . (7) 

Using  eq.  (7),  eq.  (2)  and  (3)  can  next  be  expressed  in  terms  of  the  drift 
angle  8 , 


d 6 * 

- (m+ray)U—  = Yg8+(-mU+Yr)r+Y.r+Y66  , 

(Izz«zz>£  = N68+Nrr+N&8+Nfi6  . 


(8) 

(9) 


Finally,  we  can  nondimensional ize  eq.  (5) , (8) , and  (9)  as  indicated  in  the 
Nomenclature  to  yield, 
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Y „ 1 6 ' + ( -m ' +Y  ' )r'+Y. 'r'+Y. ’6’+Y'  , (10) 

B r r <5 

N 'B'+N  1 r ' +N* 1 B ' +N  _ ' 6 ' +N 1 , (11) 

B r p 6 

*’-B'  , (12) 

where  an  external  sway  force  disturbance  Y'  and  yaw  moment  disturbance  N' 
have  been  included  for  generality. 

Equations  (10)  and  (11)  can  be  transformed  into  state-variable  form  by 
solving  eq.  (11)  for  r’  and  substituting  this  into  eq.  (10)  and  by  solving 
eq.  (10)  for  B'  and  substituting  this  into  eq.  (11).  These  can  then  be 
rearranged  to  yield. 


, . ,,dB'  _ 

-(m'+my  )-j^r  - 


dn‘ 

dt' 


dr' 

dt' 


dB' 

dt' 


22 


23 


25 


NB'(m'+my')-N.'YB' 


Nfi'(m'+my')-N&'Y6' 


N’+Y22Y'  , 

(13) 

n'+y32y'  , 

(14) 

. are  defined  as  follows: 

r 

(15) 

Y21  f26 


Yo->  ^ 


(m'+my1) 


22  27  0 ' 


32 


33 


35 


Yn  - f 


(-m'+Y r')(I  '+J  ' ) +Yf ’ Nr ' 

•r  zz  zz  L l 


Yb' (Izz*  +Jzz' )+Yr 'nB' 


Yr ' ( I ' +J  ' ) +Y* ' Ni ' 
o zz  zz  r o 


Yr' 


31  ~ 36  ~ 0 


(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 
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Y32  f 37 


(*zz  '+Jzz ' * 


(24) 


and  where, 

0=  (m'+my') (Izz’+Jz2')+Yf 'Ng'  . (25) 

Now  to  complete  the  state-variable  formulation  we  can  add  the  definition 
of  the  yaw  rate. 


d^' 

dt' 


(26) 


and  represent  the  rudder  control  system  by  a first-order  system  with  time 
constant  Tr  as  given  by  the  following  in  nondimens ional  form: 


d6' 

dt' 


(27) 


where  6C'  is  the  commanded  rudder  angle.  Equations  (12),  (13),  (14),  (26), 
and  (27)  can  now  be  rewritten  in  vector  form. 


V" 

0 1 0 0 0 

r’ 

0 f 22  f 2 3 0 f25 

B1 

se 

0 f 32  f 33  0 f 35 

n‘ 

10  -100 

_6 

00  0 0 -1/T 

V 

0 

0 

0 

r ' 

0 

^21 

^22 

B' 

+ 

0 

«'c  + 

7 31 

y 32 

V 

0 

0 

0 

_6'_ 

i/Tr 

0 

0 

(28) 


Equation  (28)  is  now  in  state-variable  form;  i.e.. 


nsxl  mxl  qxl 

x = Fs  x + Gs  u + rsw  (29) 

where  the  dimensions,  vectors  and  matrices  are  defined  by  direct  comparison 
between  eq.  (28)  and  (29) . The  subscript  s is  used  on  the  open-loop 
dynamics  matrix  Fs , control  distribution  matrix  Gs  , and  disturbance  distribu- 
tion matrix  I’s  to  signify  the  system  equations  in  contrast  to  the  augmented 
state  equations  used  below  in  the  stochastic  controller  design. 

At  this  point  it  is  useful  to  present  typical  numerical  values  for  the 

coefficients  of  the  state  equations.  The  principal  source  of  maneuvering 

coefficients  for  surface  ships  in  various  water  depths  is  the  work  of 
1 2 

Fnjino.  ’ Fujino  conducted  planar  motion  mechanism  (PMM)  and  oblique  tow 
tests  of  models  for  a Mariner  type  ship  and  for  the  tanker  Tokyo  Maru  at 
various  water  depth-to-draft  ratios,  H/T.  Table  1 gives  selected  characteris- 


tics  of  Fujino's  models.  The  Mariner  tests  were  conducted  at  Froude  numbers 
corresponding  to  7 and  12  knots  full-scale.  The  Tokyo  Maru  tests  were 
complete  only  for  a Froude  number  corresponding  to  about  10.7  knots  full- 
scale.  (Fujino  apparently  quotes  this  speed  incorrectly  as  12  knots; 
we  have  assumed  his  Froude  number  to  be  correct.)  These  data  were  used  for 
this  study.  Additional  shallow  water  data  for  a Mariner  type  ship  and  two 
Series  60  hulls  has  been  recently  presented  by  Loeser.^ 

Fujino's  data  was  transformed  using  eq.  (15)  through  (25)  and  the 
resulting  f ^ j and  y^j  coefficients  are  given  in  Tables  2,  3,  and  4. 

Table  2 presents  coefficients  for  the  Mariner  at  7 knots  full-scale  and 
values  of  H/T  equal  to  1.21,  1.50,  1.93,  2.50,  and  The  deep  water  tests 
were  conducted  at  about  H/T=30.  Table  3 presents  coefficients  for  the 
Mariner  at  12  knots  full-scale  and  values  of  H/T  equal  to  1.50,  1.93,  2.50, 
and  00 . Tests  at  H/T=1.21  were  not  feasible  at  this  speed  due  to  squat. 

Table  4 presents  coefficients  for  the  Tokyo  Maru  at  10.7  knots  full-scale 
at  values  of  H/T  equal  to  1.30,  1.50,  1.89,  2.50,  and  The  significant 
variation  in  these  coefficients  with  water  depth  and  vessel  speed  is  evident 
from  these  Tables. 


Mariner  type  ship 

Tokyo  Maru 

linear  scale  ratio,  1 

64.37 

145.0 

lenqth  between  perpendiculars,  mm 

2,500.0 

2,000.0 

breadth,  mm 

359.8 

327.6 

draft  fore,  mm 

106.5 

110.3 

aft , mm 

125.5 

110.3 

mean , mm 

116.0 

110.3 

displacement,  kg 

61.4 

58.4 

block  coefficient,  Cg 

0.5888 

0.8054 

LCB,  mm  from  amidships 

39.4  aft 

50.8  forward 

rudder  area,  mm^ 

- 

3,390.9 

propeller  diameter,  mm 

104.2 

53.8 

pitch,  mm 

108.1 

39.8 

expanded  area  ratio 

0.565 

0.619 

number  of  blades 

4 

5 

Table  1.  Characteristics  of  Fujino's  Mariner  and  Tokyo  Maru  Models 
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H/T 

1.21 

1.50 

1.93 

2.50 

OO 

f22 

-4.4390 

-2.6960 

-2.3764 

-2.1962 

-2.1939 

f23 

10.413 

7.5634 

5.6534 

4.2852 

3.3658 

^25 

-0.96338 

-1.3331 

-1.4537 

-1.5947 

-1.4609 

Y21=f26 

712.14 

822.90 

909.93 

963.51 

971.78 

>22=f27 

17.386 

14.392 

-1.8014 

-4.6587 

-4.4228 

f 32 

-0.58224 

0.012533 

0.15168 

0.18097 

0.28585 

f 33 

-1.5734 

-0.98808 

-0.80288 

-0.72058 

-0.89059 

f 35 

-0.26338 

-0.19531 

-0.23307 

-0.20998 

-0.20813 

Y31_f 36 

75.829 

30.465 

23.556 

23.294 

13.082 

Y32=f37 

-29.013 

-42.020 

-50.808 

-55.058 

-64.162 

Table  2.  Coefficients  of  Mariner  versus  H/T  at  Fn=0.0905  (7  knots  full-scale) 


H/T 

1.21 

1.50 

1.93 

2.50 

00 

f22 

★ 

-2.8841 

-2.4840 

-2.2777 

-2.3735 

f 2 3 

★ 

12.099 

6.9127 

5.2893 

3.9383 

f25 

* 

-1.1092 

-1.4273 

-1.4299 

-1.4789 

Y21=f26 

* 

844.31 

1008.1 

1042.0 

1119.6 

Y22=f27 

★ 

37.953 

-11.624 

-8.0196 

-4.2418 

f 32 

* 

-0.018982 

0.087791 

0.16516 

0.29637 

f33 

* 

-0.87677 

-0.65628 

-0.75592 

-0.80577 

f35 

* 

-0.19512 

-0.25534 

-0.19323 

-0.19941 

Y31=f36 

★ 

57.380 

75.242 

30.551 

14.603 

Y32=f 37 

* 

-29.783 

-48.715 

-52.591 

-62.167 

*not  possible  due  to  squat 

Table  3.  Coefficients  of  Mariner  versus  H/T  at  Fn=0.155  (12  knots  full-scale) 


1.30 

1.50 

1.89 

2.50 

OO 

-1.6508 

-1.7136 

-1.7657 

-1.8177 

-1.9515 

9.3157 

6.6235 

5.7359 

4.6112 

3.1591 

-0.55543 

-0.79235 

-0.88074 

-1.0416 

-1.0410 

346.69 

385.98 

477.68 

536.00 

567.13 

4.8040 

-2.2145 

-5.0043 

-5.8625 

2.3365 

0.02974 

0.13890 

0.17199 

0.23621 

0. 31507 

-1.0388 

-0.71895 

-0.52766 

-0.54560 

-0.63651 

-0.09995 

-0.12092 

-0.15607 

-0.16639 

-0.16163 

11.825 

14.230 

21.141 

21.942 

16.844 

-19.216 

-23.123 

-28.233 

-31.490 

-37.384 

Table  4.  Coefficients  of  Tokyo  Maru  versus  H/T  at  Fn=0.103  (10.7  knots  full-scale) 
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The  final  coefficient  in  eq.  (28)  requiring  numerical  definition  is  the 
nondimen sional  time  constant  of  the  rudder  control  system  Tr.  Typically, 
oceangoing  ships  have  a maximum  rudder  capability  of  ±35°  and  are  required 
to  have  a rudder  rate  capability  of  at  least  2.33°/s.  For  a first-order 
system  as  assumed  in  eq.  (27) , the  time  constant  Tr  is  the  time  at  which 
the  response  to  a step  change  in  commanded  rudder  angle  reaches  1-e  ^=0.632 
of  its  final  or  commanded  value.  With  an  initial  rudder  angle  of  zero  and  a 
commanded  35°  rudder,  it  would  require  9.5  seconds  to  reach  0.632x35=22.1° 
at  the  minimum  rate.  It  is  thus  reasonable  to  assume  a dimensional  time 
constant  of  about  10  seconds  as  done  by  Millers.^  Figure  2 illustrates  this 
situation.  The  time  constant  of  10  seconds  corresponds  to  a maximum  rudder 
rate  of  3.5  °/s.  with  a 35°  command.  Since  the  nondimensionalization  used 
here  is  speed  dependent,  a 10  second  rudder  time  constant  becomes  Tr=0.223 
for  the  Mariner  at  Fn=0.0905  or  7 knots  full-scale  and  Tr=0.382  at  Fn=0.155 
or  12  knots  full-scale.  A 10  second  rudder  time  constant  becomes  Tr=0.189 
for  the  Tokyo  Maru  at  Fn=0.103  or  10.7  knots  full-scale. 

commanded  35°  rudder 


Figure  2.  Response  of  First-order  Rudder  Model  to  '5°  Step  Command 

With  numerical  values  established  for  all  terms  in  eq.  (28)  it  is  now 
possible  to  study  the  open-loop  eigenvalues  for  the  uncontrolled  ship.  These 
were  found  using  the  OPTSYS  program  which  utilizes  Wilkinson's  QR  algorithm 
to  solve  the  open-loop  eigensystem.  The  ship  path  control  problem 
as  formulated  in  eq.  (28)  has  five  eigenvalues.  Two  eigenvalues  are  zero  and 
thus  independent  of  vessel  speed  and  water  depth.  A third  eigenvalue  is 
associated  with  the  rudder  control  system.  This  eigenvalue  equals  -1/Tr 
and  in  dimensional  form  is  also  independent  of  vessel  speed  and  water  depth. 

In  the  nondimensional  form  used  here  this  eigenvalue  appears  to  be  vessel  speed 
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dependent.  The  remaining  two  eigenvalues  vary  with  both  vessel  speed  and  water 
depth  and  are  of  primary  interest  here.  Table  5 presents  the  open-loop 
eigenvalues  for  the  Mariner  at  the  Froude  numbers  corresponding  to  7 and  12 
knots  full-scale.  The  vessel  is  course  stable  under  all  conditions.  The 
results  for  7 knots  show  an  initial  decrease  in  stability  as  shallow  water 
is  entered  (H/T=2.50)  and  then  greatly  increased  stability  in  very  shallow 
water  (H/T<1.50).  The  two  eigenvaluefe  converge  at  about  H/T=1.4  and  then 
become  a conjugate  pair  for  lower  values  of  H/T.  Table  6 presents  the  open- 
loop  eigenvalues  for  the  Tokyo  Maru  at  the  Froude  number  corresponding  to 
10.7  knots  full-scale.  The  eigenvalues  are  real  for  all  H/T.  The  tanker  is  much 
less  course  stable  than  the  Mariner  and  becomes  unstable  (eigenvalue  with 
positive  teal  part)  for  the  intermediate  range  of  water  depth-to-draft  ratio 

* 

around  H/T=2.50  and  H/T=1.89.  These  results,  of  course,  correspond  to  the 

1 ' ■ * ’ 

behavior  shown  earlier  by  Fujino  using  a different  ^approach  from  the  same  data. 


Fn=0 

.0905 

Fn=0.155 

Xl~  X2 

0 

.0 

C 

.0 

A3=-1/Tr 

-4 

.478 

*7-2 

.615 

H/T 

A4 

X5 

X4 

X5 

OO 

-.3646 

-2.720 

-.2549 

-2.924 

2.50 

- . 3095 

-2.607 

-.3116 

-2.722 

1.93 

-.3745 

-2.805 

-.3693 

-2.771 

1.50 

-.9343 

-2.750 

-.9986 

-2.762 

1.21 

-3.006  ± 

2 . 002  j 

not  feasible 

Table  5.  Open-loop  Eigenvalues  for  Mariner 


Xl=  X2 

0.0 

A3=-1/Tr 

-5.281 

H/T 

a4 

X5 

OO 

-.0992 

-2.489 

2.50 

.0405 

-2.404 

1.89 

.0237 

-2.317 

1.50 

-.1358 

-2.297 

1.30 

-.7360 

-1.954 

Table  6.  Open-loop  Eigenvalues  for  Tokyo  Maru  at  Fn=0.103 


The  open-loop  eigenvectors  for  the  uncontrolled  ship  are  also  of  interest. 
The  repeated,  zero  eigenvalue  has  only  one  linearly  independent  eigenvector 
[0,0, 0,1,0]  ; i.e.,  just  n'-  This  can  be  viewed  as  the  steady-motion 
condition  where  the  ship  has  constant  lateral  offset;  i.e.  r ' = 6 '=ij/ ' = 6 ’=0 
but  n V0.  A second,  generalized  eigenvector^  can  be  associated  with  the 
zero  eigenvalues.  This  eigenvector  can  be  shown  to  be  [a, 0,0, 1,0]  where  a 
is  any  nonzero  quantity.  This  can  be  viewed  as  the  steady-motion  condition 
where  the  ship  has  a constant  heading  away  from  the  desired  track;  i.e., 
r'=8'=6'=0  but  i|>V0  and  n VO.  The  remaining  three  eigenvectors  associated 
with  X3,  X4,  and  X5  are  vessel  speed  and  water  depth  dependent.  These 
eigenvectors  are  shown  for  the  Mariner  at  four  water  depth  and  vessel  speed 
conditions  in  Table  7 and  for  the  Tokyo  Maru  at  Fn=0.103  in  Table  8.  The 
eigenvector  associated  with  X3  is  the  only  one  containing  a 6'  component 
as  expected.  This  eigenvector  has  some  water  depth  dependence  but  strong 
vessel  speed  dependence  as  shown  in  Table  7 due  primarily  to  the  improved 
rudder  effectiveness  with  vessel  speed.  The  eigenvectors  associated  with 
X4  and  X5  show  much  smaller  vessel  speed  dependence  but  a much  larger  water 
depth  dependence  for  both  types  of  ships.  The  eigenvector  associated  with 
the  dominant  eigenvalue  X4  shows  a major  change  in  character  as  the  ships 
become  extremely  course  stable  in  very  shallow  water. 
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eigenvector  associated  with  A ,=  -1/Tr 


Fn 

0.0905 

0.155 

0.0905 

0.155 

H/T 

CO 

OO 

1.50 

1.50 

ip  ' component 

-0.1179 

0.3300 

-0.1021 

-0.3138 

r ' component 

0.5278 

-0.8629 

0.4571 

0 . 8205 

3 ' component 

0.0067 

0.1782 

0.0477 

0.0597 

ri 1 component 

0.0278 

-0.0581 

0.0335 

0.1429 

6 ' component 

0.8406 

0.3338 

0.8816 

0.4521 

eigenvector  associated  with  A 

4 

Fn 

0.0905 

0.155 

0.0905 

0.155 

H/T 

CO 

OO 

1.50 

1.50 

0 ' component 

-0.2890 

-0.2813 

-0.5257 

-0.5455 

r component 

0.1054 

0.0556 

0.4911 

0.5447 

3 component 

0.0573 

0.0299 

0.1144 

0.0849 

h component 

0.9498 

0.9738 

0.6851 

0.6313 

6 ' component 

0.0000 

0.0000 

0.0000 

0.0000 

eigenvector  associated  with  A 

5 

Fn 

0.0905 

0.155 

0.0905 

0.155 

H/T 

oo 

OO 

1.50 

1.50 

tp ' component 

-0.3405 

-0.3201 

-0.3392 

-0.3377 

r ' component 

0.9262 

0.9361 

0.9329 

0.9328 

3 ' component 

-0.1447 

-0.1309 

-0.0066 

0.0094 

n ' component 

0.0720 

0.0647 

0.1210 

0.1256 

6 1 component 

0.0000 

0.0000 

0.0000 

0.0000 

Table  7.  Open-loop  Eigenvectors  for  Mariner  at  Four  Conditions 
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eigenvector  associated  with  A3  = -1/Tr 


H/T 

CO 

2.50 

1.89 

1.50 

1.30 

component 

-0.0541 

-0.0496 

-0.0387 

-0.0341 

-0.0178 

r ' 

component 

0.2859 

0.2621 

0.2046 

0.1801 

0.0938 

8' 

component 

0.0139 

0.0208 

0.0247 

0.0206 

0.0228 

n ' 

component 

0.0129 

0.0133 

0.0120 

0.0104 

0.0077 

6 1 

component 

0.9566 

0.9635 

0.9777 

0.9828 

0.9951 

eigenvector  associated  with 

a4 

H/T 

OO 

2.50 

1.89 

1.50 

1.30 

** 

component 

-0.0933 

0.0412 

0.0238 

-0.1304 

-0.5220 

r ' 

component 

0.0093 

0.0017 

0.0006 

0.0177 

0.3842 

8' 

component 

0.0054 

0.0007 

0.0002 

0.0042 

0.0377 

n' 

component 

0.9958 

0.9991 

0.9997 

0.9913 

0.7606 

6' 

component 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

eigenvector  associated  with  Ag 


H/T 

OO 

2.50 

1.89 

1.50 

1.30 

r 

component 

-0.3669 

-0.3792 

-0.3913 

-0.3942 

-0.4450 

r' 

component 

0.9132 

0.9115 

0.9067 

0.9053 

0.8693 

8' 

component 

-0.1553 

-0.1159 

-0.0871 

-0.0797 

-0.0283 

n' 

component 

0.0850 

0.1095 

0.1313 

0.1369 

0.2133 

6 ' 

component 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Table  8.  Open-loop  Eigenvectors  for  Tokyo  Maru  at  Fn=0.103 


3.  Optimal  Stochastic  Path  Controller  Design 


This  section  will  begin  with  a brief  introduction  to  the  optimal  control 
of  stochastic  linear  systems.  Readers  already  familiar  with  these  ideas  can 
continue  to  the  second  subsection. 
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3.1  Introduction  to  Optimal  Stochastic  Control  ' ' 

Many  physical  control  problems  can  be  realistically  represented  using 
stochastic  disturbances  and  measurement  noise.  This  is  particularly  true  where 
disturbances  are  of  a random  nature  or  where  they  cannot  be  defined  specifically 
in  advance.  The  response  of  a system  subjected  to  stochastic  disturbances  and 
measurement  noise  can  also  be  treated  as  a stochastic  quantity.  The  engineering 
approach  can  be  to  model  stochastic  physical  systems  as  Gauss-Markov  pro- 
cesses. These  can  be  represented  by  the  state  vector  of  a linear  dynamical 
system  forced  by  a gaussian  purely-random  process  where  the  initial  state 
vector  is  also  gaussian  or  normally  distributed.  Thus  we  can  represent  the 
system  by  eq.  (29)  which  we  repeat  again  here  for  reference;  i.e., 

nsxl  mxl  qxl 

x = Fsx  + Gs  u + rsw  . (29) 

The  system  open-loop  dynamics  matrix  Fs  , the  system  control  distribution 
matrix  Gs  , and  the  disturbance  distribution  matrix  ?s  are  assumed  constant 
here.  The  condition  of  this  system  is  completely  represented  by  the  mean 
value  vector  x and  the  covariance  matrix  X for  the  ns  differentiated  or 
state  variables;  i.e., 

x(t)  = E [x (t) ] with  x(tQ)  = xQ  , 

X (t)  = E [ (x (t) -x  (t) ) (x  (t) -x (t) )T]  with  X(tQ)=X0  . 

where  E ( . . . ] is  the  expected  value  or  ensemble  average  over  the  many 
possible  observations  at  time  t.  The  m non-dif ferentiated  variables  in 
u are  the  control  variables.  The  q variables  in  w are  the  process 
disturbances  which  are  gaussian  purely-random  processes  or  white  noise. 

White  noise  is  an  idealized,  very-jittery  process  which  can  be  viewed  as  the 
limit  of  a sequence  of  impulses  with  random  magnitude  and  random  time  of 
occurrence.  The  impulses  average  zero  over  time  but  have  an  average  square 
magnitude  given  by  o(t)  squared.  We  thus  have, 


(30) 

(31) 
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E [ w ( t ) ] = 0 , (32) 

E[w(t)w(x)T]  = Q(t)6(t-T)  , (33) 

where  Q is  the  power  spectral  density  matrix  which  can  be  related  to  o(t) 
and  S(t-i)  is  the  Dirac  delta  function.  It  is  also  assumed  that  there  is 
no  correlation  between  the  process  disturbance  and  the  initial  condition  of 
the  system;  i.e. 

E(w(t) (x(t0)-x0)T)  = 0 . (34) 

The  control  problem  is  to  develop  the  optimal  state-variable  feedback 
control , 


u = Cx 


where  C is  the  control  feedback  gain  matrix.  In  general,  not  all  of  the 
needed  states  in  x are  readily  measured.  Further,  it  is  not  necessary  to 
measure  all  the  states  if  it  is  possible  to  estimate  the  remaining  states  from 
those  which  are  most  easily  measured.  In  the  stochastic  case  we  may  have  p 
available  measurements  represented  by, 

pxl 

1 = Hs*  + v ' (36) 

where  Hs  is  the  measurement  distribution  matrix  (assumed  constant  here) 
and  v is  a vector  of  white  measurement  noise  with  statistical  properties, 

E[v(t) 1=0  , (37) 

E[v(t)v(T)T]  = R (t ) 6 ( t-T ) , (38) 

E ( v (t ) w ( t) T]  = E[v(t) (x(t0)-^0)T]  = 0 . (39) 

The  matrix  R is  the  power  spectral  density  of  the  measurement  noise. 

Equation  (39)  states  that  there  is  no  correlation  between  the  measurement 
noise  and  the  process  disturbance  or  the  initial  state  of  system.  The  elements 
of  z may  be  measurements  of  specific  states  or  linear  combinations  of  the 
states. 


In  modeling  physical  systems,  it  is  not  always  realistic  to  assume  that 
the  process  disturbance  or  measurement  noise  is  white  noise.  If  the  process 
disturbance  is  a random  quantity  which  changes  very  rapidly  compared  to  the 
time  response  of  the  system,  it  is  reasonable  to  assume  the  process  disturbance 
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is  white  noise  as  in  Eq.  (29)  . However,  if  the  process  disturbance  is  a 
random  quantity  which  changes  very  slowly  compared  to  the  time  response  of 
the  system  (perhaps  a tidal  current  effect  on  a passing  ship) , it  is  reasonable 
to  assume  the  process  disturbance  to  be  a random  bias  or  constant.  This  can 
be  incorporated  into  the  above  treatment  by  defining  an  additional  state 
variable  or  variables  such  that, 

xn+i  = 0 , (40) 

and  xn+i (tQ)  is  random.  This  state  variable  is  included  in  an  augmented 
state  vector  of  length  n+1  with  the  n+1  component  from  the  disturbance  equal 
to  zero. 

If  the  process  disturbance  is  a random  quantity  which  changes  on  about 
the  same  time  scale  as  the  response  of  the  system,  it  must  be  modeled  as 
something  between  white  noise  and  a random  bias.  This  is  accomplished  by  the 
use  of  various  shaping  filters  and  again  augmenting  the  state  vector.  The 
simplest  shaping  filter  produces  an  exponentially  correlated  disturbance^ 
by  driving  a first-order  system  by  white  noise.  A new  state  variable  is 
defined  as  follows: 

T *n+l  + xn+l  = w ' <41> 

where  T is  the  disturbance  correlation  time  and  w is  white  noise.  The 

state  vector  can  then  be  augmented  to  an  n+1  vector  and  the  total  system 

is  still  disturbed  by  white  noise  as  in  eq.  (29) . The  shaping  filter  processes 

the  white  noise  to  produce  a new  disturbance  xn+i  which  is  random  but  with 

a characteristic  time  constant  t of  about  the  same  order  as  the  response 

time  of  the  system.  Other  higher-order  shaping  filters  are  possible  to  model 

3 5 18 

more  complex  disturbances.  ' ' 

If  shaping  filters  are  used  to  mode)  some  ne-ns  of  the  disturbances,  the 
augmented  state  vector  is  then  of  length  ne  and  the  stochastic  controller 
design  is  performed  using  a set  of  augmented  or  estimator-design  state 
equations  combining  eq.  (29)  and  ne-ns  equations  like  eq.  (41);  i.e., 

nexl  mxl  qxl 

x’  = Fex’  + Ge  u + re  w (42) 

where  the  disturbance  w is  composed  of  only  white  noise  with  a power  spectral 
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density  given  by  Q.  To  complete  the  formulation,  the  measurement  eq.  (36)  can 
now  be  expressed  in  terms  of  the  augmented  state  vector  to  be. 


z = (Hc 


0)  x ' + v = Hex'  + v 


(43) 


The  Separation  Theorem  states  that  the  optimal  way  to  control  the  system 
eq.  (29)  which  has  been  modeled  by  eq.  (42)  using  the  information  available 
in  the  noisy  measurements  eq.  (43)  is  to  design  controller  gains  neglecting 
w and  v and  thus  assuming  perfect  knowledge  of  x'  . This  control  would  be 
given  by. 


u = Cx‘ 


(44) 


The  noisy  measurements  z_  can  then  be  utilized  in  an  optimal  stochastic 
observer  (state  estimator)  or  Kalman-Bucy  filter  to  produce  a maximum- 
likelihood  estimate  of  the  augmented  state  vector  x which  includes  estimates 
of  the  actual  states  x plus  estimates  of  any  augmented  states  created  in  the 
modeling.  This  estimated  state  vector  is  available  to  be  used  in  eq.  (44) 
so  the  optimal  control  will  then  be  given  by, 


u = Cx 


(45) 


Thus,  the  controller  design  can  be  completely  separated  from  the  processing 
of  the  noisy  measurements  to  produce  the  best  estimate  of  the  current  state 
of  the  system.  An  overall  schematic  of  such  a stochastic  control  system  is 
shown  in  Figure  3.  The  measurements  from  the  sensors  are  used  in  the 

optimal  stochastic  observer  to  produce  an  estimate  of  the  states  k which 
are  then  used  in  the  optimal  controller  to  produce  the  control  signals  ij 
given  to  the  actuators. 

An  optimal  control  can  be  defined  in  many  ways.  The  most  common  when 
we  want  to  control  x near  zero  using  reasonable  values  of  control  u is 
to  use  the  control  which  minimizes  a linear  quadratic  cost  function, 


= 2 j (x'TAx'  + uTBu! 


Bu)  dt 


(46) 


where  the  A and  B matrices  are  initially  established  by  the  designer  to 
reflect  the  relative  weighting  of  errors  in  the  various  states  and  the  use 
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Figure  3.  Overall  Schematic  of  Optimal  Stochastic  Control  System 

of  the  various  controls.  The  A and  B matrices  are  usually  diagonal  with 
at  least  some  nonzero  diagonal  elements  selected  to  be, 


Aii  x .2 

xoi 


and  Bjj  = 


where  uQj  is  an  acceptable  amount  of  control  j to  be  used  when  the  state 
i deviates  x0^  from  zero.  It  is  usually  necessary  to  modify  the  weighting 
matrices  A and  B and  iterate  on  the  design  based  on  the  evaluated  response 
of  the  controlled  system.  The  calculus  of  variations  ' can  be  utilized  to 
show  that  the  control  which  minimizes  eq.  (46)  for  a stationary  system  is 
given  by, 

C = , (48) 

6 00 

where  S is  the  steady-state  solution  of  a matrix  Riccati  equation  which  is 

OO 

independent  of  the  process  disturbance  w and  Q and  the  measurement  noise 
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v and  R ; i . e . , 

S = -SFe-FeTS+SGeB_1GeTS-A  . (49) 

The  open-loop  dynamics  matrix  Fg  and  the  control  distribution  matrix  Ge 

are  for  the  augmented  state  equations  (42)  actually  used  in  the  controller 

design.  If  shaping  filters  are  not  used  in  the  modeling  of  the  process 

disturbances,  eq.  (29)  and  (42)  are  identical  so  ne=ns , Fe=Fs,  Ge=Gs , and 

re=rs  • An  efficient  way  to  obtain  the  steady-state  solution  to  eq.  (49), 

, is  to  utilize  the  technique  of  eigenvector  decomposition  first  proposed 
19  20 

by  MacFarlane  and  Potter  . This  method  was  developed  into  a practical 
design  computer  program  by  Bryson  and  Hall'*'0  in  Hall's  OPTSYS  program  which 
uses  the  QR  algorithm  to  solve  the  eigensystem.  User's  documentation  is 
available*'*'  for  the  Michigan  version  of  the  OPTSYS  program  as  used  in  this 
study. 

The  second  half  of  the  design  problem  is  to  develop  the  optimal  stochastic 

17 

observer.  Again,  the  calculus  of  variations  can  be  utilized  to  show  that 
the  maximum-likelihood  estimate  of  the  state  is  produced  by  the  filter  given 
by, 

k = Fe*  + GeH.  + K(z-Hex)  ; x(tD)  = xQ  , (50) 

where  the  filter  gain  matrix  K for  a stationary  system  is  given  by, 

K = P^Hg^R-1  , (51) 

where  P^  is  the  steady-state  solution  of  the  matrix  Riccati  equation, 

P = FeP  + PFeT  + reQreT  - PHeTR-1HeP  . (52) 

The  matrix  P is  the  covariance  of  the  error  of  the  estimate  of  the  state; 
i .e.  , 

P(t)  = [E  (x  (t) -x*(t) ) (x(t)-x'(t)  )T]  . (53) 

In  eq.  (50),  it  can  be  seen  that  the  estimate  x is  assumed  to  follow  the 
same  dynamics  as  x'  (excluding  IV)  and  that  whenever  the  measurement  which 
would  result  from  k (excluding  v)  deviates  from  the  actual  measurement  z 
a correction  is  introduced  to  drive  k closer  to  x'  . In  eq.  (52) , it  can 
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be  seen  that,  as  expected,  large  process  disturbances  (high  Q)  and  large 
measurement  noise  (high  R)  will  increase  P,  the  error  in  the  estimate  produced 
by  the  stochastic  observer.  An  efficient  way  to  obtain  the  steady-state  solu- 
tion to  eq.  (52)  , P , is  again  by  eigenvector  decomposition  as  implemented 
in  the  OPTSYS  program  used  in  this  study. 

Summarizing,  if  a physical  system,  its  disturbances,  and  measurement 
noise  are  modeled  as  described  above,  a design  program  such  as  the  OPTSYS 
program  can  be  used  to  produce  the  optimal  control  gains  C and  the  optimal 
filter  gai  s K.  The  OPTSYS  program  can  also  produce  the  Root  Mean  Square 
(RMS)  response  of  the  system  to  the  design  disturbances  as  represented  in 
the  power  spectral  densities  Q and  R.  The  RMS  response  is,  however,  not  that 
meaningful  to  many  engineers.  Further  if  specific  physical  disturbances  were 
modeled  in  the  design  process  using  shaping  filters,  the  designer  may  also 
want  to  know  the  response  of  the  controlled  system  to  the  specific  disturbances . 
Thus,  there  is  often  a need  to  simulate  the  response  of  the  optimal  stochastic 
control  system  to  specific  process  disturbances  and  initial  condition  errors 
while  subject  to  the  measurement  noise.  The  SHIPSIM/OPTSIM  simulation 
program^  which  we  have  developed  as  a complement  to  the  OPTSYS  program 
allows  the  simulation  of  these  systems  with  a minimum  of  effort  and  computer 
programming . 

3. 2 Process  Disturbance  Modeling 

A ship  moving  along  a prescribed  path  can  be  subjected  to  many  types  of 

short  term  disturbances.  Waves  provide  a relatively  high  frequency  disturbance. 

Changes  in  current,  banks,  and  bottom  can  also  provide  significant  disturbances. 

For  this  study,  we  selected  as  the  design  disturbance  the  yaw  moment  and 

sway  force  due  to  a passing  ship  which  constitute  one  of  the  strongest 

disturbances  a ship  can  experience  in  maneuvering  situations.  Data  given  by 
21 

Newton  was  utilized  to  establish  a representative  magnitude  and  time  history 

for  this  disturbance.  His  particular  data  were  for  a Mariner  moving  in  deep 
22  2 3 

water.  Yung  and  Abkowitz,  Ashe,  and  Fortson  show  a significant  increase 
in  the  disturbance  yaw  moment  and  sway  force  due  to  a passing  ship  as  water 
depth  decreases  so  Newton’s  data  is  clearly  lower  than  would  be  expected  in 
shallow  water.  For  the  purposes  of  this  study,  we  wanted  only  a representative 
disturbance  so  the  yaw  moment  and  sway  force  design  disturbances  shown  in 


Fig.  4 were  chosen  and  assumed  to  be  independent  of  ship  type,  and  water  depth, 
and  vessel  speed. 


Figure  4.  Design  Yaw  Moment  and  Sway  Force  Disturbances  due  to 

Passing  Ship 


Disturbances  as  assumed  in  Fig.  4 can  be  reasonably  modeled  as  exponen- 
tially correlated  disturbances  as  produced  by  first-order  shaping  filters 
driven  by  white  noise.  Thus,  two  new  states  were  defined  by. 


(54) 

(55) 


where  the  correlation  times  were  taken  as  TN=Ty=l  or  dimensionally  the 
time  it  takes  the  ship  to  travel  its  own  length.  The  power  spectral  densities 
for  the  white  noise  wN  and  wy  were  defined  based  on  the  Root  Mean  Square 

(RMS)  values  of  N'  and  Y',  respectively,  during  the  t'  period  (-2, +1.5) 

- -4  -4 

shown  on  Fig.  4;  i.e.,  N' =0.880x10  and  Y' =2. 118x10  . The  white  noise  power 

spectral  densities  are  then  given  by3. 


q 


11 


2 (N ' ) 


-8 


1.548x10 


(56) 
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for  wN  and, 

q22  “ 2(?')2Ty=8.970x10~8  , (57) 

for  wy.  The  design  disturbances  are  fully  defined  by  q^,  q22»  qi2=q21=°  ' 
and  eq.  (54)  and  (55).  Simulation  of  the  controlled  ship  subjected  to  the 
specific  disturbances  shown  in  Fig.  4 have  confirmed  the  validity  of  this 
modeling  approach  as  will  be  discussed  below. 


With  two  shaping  filters  used  to  model  the  yaw  moment  and  sway  force 
disturbances,  the  augmented  system  eq.  (42)  used  in  the  stochastic  control 
system  design  combine  eq  (28),  (54),  and  (55)  to  give, 
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(58) 


The  dimensions,  vectors,  and  matrices  of  eq.  (42)  are  defined  by  direct 
comparison  with  eq.  (58).  Notice  that  the  upper  right  (5x2)  partition  of 
Fe  in  eq.  (58)  is  Ts  in  eq.  (28);  i.e.  f26=Y21'  Y36=Y31'  etc‘ 


3 . 3 Selection  of  Measurements 

The  measurements  eq.  (43)  and  associated  measurement  noise  power  spectral 
density  matrix  R were  established  next.  In  this  study,  we  considered  the 
use  of  a ip'  measurement  from  a compass  or  gyro  compass,  an  r1  measurement 
from  a rate  gyro,  a f?'=-v'  measurement  from  a doppler  sonar,  and  an  n' 
measurement  from  radar  or  DECCA.  Rudder  angle  could  be  easily  obtained  but 
was  not  considered  here.  The  first  consideration  to  answer  in  selecting 
measurements  is  whether  or  not  all  the  states  of  the  path  control  problem  are 
observable  with  a given  measurement  or  set  of  measurements  z.  If  a state 
is  not  observable  with  a given  , it  is  impossible  to  design  a filter  which 
can  produce  the  augmented  state  vector  estimate  x from  the  measurement 
vector  z_. 

One  way  to  establish  the  observability  of  eq.  (42)  and  eq.  (43)  is  through 

16 

the  use  of  a Jordan “form  transformation  . If  the  open- loop  eigenvectors  of 
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Fg  are  grouped  together  as  columns  of  a transformation  matrix  T,  the  problem 


can  be 
£'  =T_1x 

transformed  into  its  equivalent  Jordan-form  using  the  transformation 

: . Using  this  transformation  eq.  (42)  and  (43)  become, 

£’ 

= A£  + G' 

1 u + r • w , 

(59) 

and , 

z 

II 

X 

+ 

< 

$ 

(60) 

where. 

A 

= T_1FeT 

t 

(61) 

G' 

= T_1Ge 

t 

(62) 

T1 

= T_1re 

r 

(63) 

H' 

II 

X 

rt> 

• 

In  this  form,  the  various  modes  of  response  of  the  system  are  completely 
decoupled  provided  the  eigenvalues  of  Fg  are  distinct.  In  that  case,  matrix 
A is  a diagonal  matrix  with  the  open-loop  eigenvalues  of  Fe  down  the 
diagonal.  The  structure  of  H'  can  be  studied  to  quickly  establish  the 
observability  of  each  mode  of  the  system  with  each  measurement  in  z_.  Similarly, 
the  structure  of  G'  can  be  studied  to  establish  the  controllability  of  each 
mode  of  the  system  with  each  component  of  the  control  vector  u. 

The  ship  path  control  problem  as  formulated  in  eq.  (58)  has  seven  eigen- 
values, two  of  which  are  zero,  but  only  six  independent  eigenvectors.  The 

Jordan-form  transformation  can  be  performed  using  these  six  eigenvectors  plus 

T 

an  additional,  generalized  eigenvector  (i  0 0 1 0 0 0 ] , a>0  as  discussed  in 
Section  2.  Four  of  the  eigenvectors  are  identical  to  those  presented  in 
Section  2 with  two  zeros  added  to  make  them  of  dimension  (7x1).  Two  additional 
eigenvectors  are  associated  with  the  shaping  filters.  The  essential  elements 
of  the  resulting  Jordan -form  structure  are  shown  in  Fig.  5 where  C 2'  is 
associated  with  the  generalized  eigenvector  for  the  repeated  eigenvalue 
Ai=A2=0  . The  OPTSYS  program  was  utilized  to  perform  the  Jordan-form  trans- 
formation and  to  establish  numerical  values  for  matrices  G' , H',  and  T'.  If 
Hji'  is  zero  it  shows  directly  that  mode  is  not  observable  with  measure- 

ment j,  etc.  Likewise,  if  G31'  is  zero  it  directly  shows  that  mode  £3' 
cannot  be  controlled  with  the  single  control  6^'  , etc. 
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Figure  5.  Simplified  Jordan-Form  Structure  of  Path  Control  Problem 

The  OPTSYS  program  also  presents  the  elements  of  H'  in  a normalized 
manner  based  on  the  orthogonality  between  the  eigenvector  of  T associated 
with  mode  i and  the  row  of  H associated  with  measurement  j.  A result 
cos6^=l  indicates  that  mode  i is  fully  observable  with  measurement  j ; a 
result  cos?  ^=0  indicates  that  mode  i cannot  be  observed  (or  estimated) 
based  on  measurement  j.  The  numerical  value  between  0 and  1 thus  indicates 
quantitatively  how  observable  each  mode  is  with  any  particular  measurement. 
The  observability  results  for  the  Tokyo  Maru  at  H/T=1.89  and  Fn=0.103 
are  shown  in  Table  9.  This  shows  that  mode  is  only  observable  with  a 

measurement  of  n ' so  as  a minimum  n ' must  be  measured.  A measurement  of 
0*  also  gives  good  information  on  modes  £ ' and  and  will  allow 

estimation  of  the  sway  force  disturbance  represented  in  mode  £7'.  The 
observability  of  the  other  modes  with  n'  is  much  weaker  so  at  least  one  more 
measurement  is  probably  necessary.  Adding  anr'  measurement  gives  good  obser- 
vability of  modes  ' » £5'  and  the  yaw  moment  disturbance  represented  in 
mode  Thus,  measurement  of  r'  and  if  would  probably  be  sufficient.  In 

our  study,  we  have  used  measurements  of  , r',  and  n'  as  the  base  set  of 


mode 

*1’ 

0.0000 

0.0000 

0.0000 

1.0000 

mode 

f (O0*0 

0.0000 

0.0000 

1.0000 

mode 

h' 

0.0238 

0.0005 

0.0002 

0.9997 

mode 

W 

0.3913 

0.9067 

0.0872 

0.1313 

mode 

*5* 

0.0388 

0.2046 

0.0247 

0.0120 

mode 

^(N') 

0.5881 

0.5881 

0.5542 

0.0339 

mode 

cytY') 

0.5487 

0.5487 

0.0773 

0.6260 

Table  9.  Observability  cos(3.  . of  Tokyo  Maru  Modes  with 
Various  Measurements 

measurements  but  Table  9 shows  that  the  measurement  could  probably  be 

omitted  to  simplify  the  control  system  without  serious  loss  of  controller 
performance.  Table  9 also  shows  that  the  use  of  a doppler  sonar  to  produce 
a drift  angle  measurement  f? ' could  probably  not  be  justified. 

Using  the  measurements  of  ip',  r',  and  g ' eq.  (43)  becomes. 


z 


1 0 0 0 0 0 0 

0 1 0 0 0 0 o 

0 0 o 1 o 0 0 


x'  + V 


(65) 


The  measurement  noise  power  spectral  density  matrix  R for  the  vector  v can 
be  estimated  for  each  measurement  by  assuming  an  RMS  noise  value  Oj  and  a 
correlation  time  Tj  (much  faster  than  time  constants  of  the  ship) . The 
white  noise  power  spectral  density  component  for  measurement  is  then  given  by, 

r^.  j = 2 (Oj)2x j . (66) 

The  numerical  values  for  the  noise  power  spectral  densities  of  the  four 
possible  measurements  used  in  this  study  are  summarized  in  Table  10.  In  the 
nondimensional  form  shown  in  Table  10  the  resulting  power  spectral  densities 

r..  are  different  for  each  ship  and  speed.  Data  Oj  and  Tj  for  the 

J 9 

measurements  of  t|/',  r',  and  g'  are  from  Millers.  Data  for  the  doppler 
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sonar  are  estimated  based  on  information  for  the  Furuno  MF-100DS.  Using 
measurements  of  ip',  r',  and  n'  as  in  eq.  (65),  the  matrix  R is  diagonal 


J 


with  associated  diagonal  elements  taken  from  Table  10. 
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r . . 

D D 

measurement 

source 

0 . 

I . 

D 

Mariner 

Mariner 

Tokyo  Maru 

Fn=0.0905 

Fn=0.155 

Fn=0 . 103 

gyro  compass 

0.2° 

0.  Is. 

5.441x10“^ 

9. 320xl0”8 

4.612xl0~8 

r 1 

rate  gyro 

'J.01°/s 

0.1s. 

2.728xl0”7 

1. 593xl0_7 

3.218xl0~7 

6* 

doppler 

sonar 

0 . 05m/s 

0.1s. 

8.629xl0“7 

5.007xl0"7 

3 . 1 39xl0~  7 

n' 

i adar 

10m 

0.1s. 

1.725xlO-5 

2.955xl0_S 

4.502xl0"6 

Table  10.  Measurement  Noise  Power  Spectral  Densities 


The  ship  path  control  problem  as  formulated  in  eq.  (58)  is  known  to  be 
controllable  with  the  rudder.  The  G'  matrix  established  for  this  problem 
using  the  OPTSYS  program  verifies  this  by  having  nonzero  elements  G^',  G3i ' ' 
, and  1 which  indicate  that  all  the  real  states  of  the  system  can  be 
controlled  using  the  rudder  command  6 '.  As  shown  in  Fig.  5,  elements  G, . ' 

C bl 

and  Gyi'  are  zero  indicating,  again  as  expected,  that  the  rudder  command  6c’ 
has  no  influence  on  the  yaw  moment  and  sway  force  disturbances  represented 
in  mode  and  mode  Cy'f  respectively. 


3 . 4 Design  Results  for  the  Tokyo  Maru  at  H/T=1.89 

The  example  results  to  be  presented  in  this  section  are  for  the  Tokyo 
Maru  at  Fn=0.103  and  H/T=1.89.  As  will  be  shown  in  the  next  section,  if 
the  stochastic  path  control  system  is  to  be  designed  for  ship  characteristics 
at  a single  depth  it  is  best  to  choose  one  of  the  depths  at  which  the  vessel 
is  least  course  stable.  The  Tokyo  Maru  is  course  unstable  in  the  approximate 
range  1 . 75<h/T<2 . 75  so  a design  depth  of  H/T=1.89  was  selected  for  the  base 
case  to  be  presented  here.  The  system  is  fully  defined  by  eq.  (58) , eq. 

(65),  and  the  various  data  presented  above.  The  only  remaining  requirement 
is  to  establish  the  state  weighting  matrix  A and  the  control  weighting 
matrix  B used  in  eq.  (46)  to  define  the  basis  for  determining  the  optimal 
control  gains  C. 


For  the  ship  path  control  problem,  we  want  to  control 

zero  without  using  excessive  amounts  of  control  6 ' 

c 


n ' and  ip  ' 


near 


and  thus  also  6'. 


For 
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this  study,  we  assumed  that  we  would  be  willing  to  commit  5°  rudder  usage 
for  path  control  when  the  ship  heading  deviates  5°  or  the  ship  reaches  a 
lateral  offset  ri0  of  10.4  3m  for  the  Tokyo  Mar u or  5.79m  for  the  Mariner. 

IThis  lateral  offset  is  exactly  a quarter  beam  (B/4)  for  the  Mariner  and 

slightly  less  than  a quarter  beam  for  the  Tokyo  Maru . A 5°  heading  deviation 
for  either  ship  will  cause  the  ends  of  the  ship  to  be  about  the  same 
distance  from  the  desired  track  as  the  assumed  nQ  value  (12.7m  for 
Tokyo  Maru ; 7.0m  for  Mariner) . Using  eq.  (47),  the  A and  B matrices  were 
taken  as  zero  except  for  the  following  elements : 


> = (5/57. 3)“2  = 131.332  for  ip ' , 

(67) 

= (10.43/290)  =772.463  for  n. 

(68) 

55  = (5/57.3)  = 131.332  for  6'  , 

(69) 

i = B,  =(5/57. 3)~2=131. 332  for  6 1 . 

11  c 

(70) 

These  A and  B matrices  were  used  for  all  Tokyo  Maru  and  Mariner  designs 
presented  here.  These  state  deviation  and  control  weighting  terms  can  be 
modified  at  the  designer's  discretion  based  on  the  performance  of  the 
controlled  system.  Increasing  the  value  of  a term  in  the  A matrix  will 
produce  a system  design  which  will  more  closely  control  the  associated 
state  at  the  expense  of  increased  deviation  of  the  other  states  and  the  use 
of  more  control.  Increasing  the  value  of  B will  reduce  the  amount  of  commanded 
rudder  (and  also  6')  at  the  expense  of  greater  deviations  in  the  other 
states.  For  the  purposes  of  the  study  presented  here,  the  A and  B matrices 
were  held  constant  (following  preliminary  investigations)  to  provide 
meaningful  comparisons  among  competing  designs. 

Using  one  control  6c'  to  control  five  states  simultaneously,  it  is 
difficult  to  directly  compare  the  RMS  response  of  the  ship  to  the  design 
disturbances  when  using  various  control  systems.  We  have  therefore  used 
an  additional  control  measure-of-merit  which  we  call  the  RMS  cost  J.  This 
quantity  is  based  on  eq.  (46)  and  reflects  the  weighting  implied  by  eq. 

(67)  through  eq.  (70).  The  RMS  cost  is  defined  by, 

j = ^,2+(A44n'2+A555'2+B115c'2)/A11  , (71) 


where  (.7.)  indicates  the  RMS  value  of  the  particular  quantity.  This  RMS 
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cost  J is  just  twice  the  ~ost  J defined  by  eq.  (46)  divided  by  A^. 

For  the  numerical  values  used  here  this  becomes, 

J = ij),2+5.8818n'2+6,2+6c'2  . (72) 

The  cost  J or  J provides  the  most  meaningful  single  quantity  with 
which  to  compare  the  RMS  response  of  a ship  when  controlled  by  two  different 
control  systems  which  were  developed  using  the  same  state  and  control  weighting 
matrices . 

Using  eq.  (58) , eq.  (65) , and  the  A and  B weighting  matrices  defined 
above,  the  OPTSYS  program  was  used  to  determine  the  optimal  control  gains 
C and  the  optimal  filter  gains  K for  the  Tokyo  Maru  at  Fn=0.103  and  H/T= 

1.89.  The  resulting  gains  are  as  follows: 


” 5.6470~ 

~ 2.6321 

0.9164 

0.0015" 

2.6848 

6.3936  12.2380 

-0.0724 

6.3147 

1.4384 

3.2622 

-1.2160 

2.4252 

; K = 

0.1447  -1.0125 

1.5367 

-0.8079 

0.0000 

0.0000 

0.0000 

675.14 

0.0263 

0.1673 

0.0130 

-50.154 

-0.0287  -0.1710 

0.0319 

Notice  that  the  estimates  of  the  yaw  force  disturbance  N'  and  the  sway 
force  disturbance  Y'  produced  by  the  filter  are  included  in  the  augmented 
state  feedback;  i.e.,  C.^0,  C,-,^0.  Also  notice  that  the  Kalman-Bucy  filter 

lb  1 / 

does  nothing  (fifth  row  of  K is  zero)  to  improve  the  existing  estimate  of  the 
rudder  angle  6'.  This  is  because  6'  is  not  measured  and  its  response  depends 
only  on  the  current  6'  and  the  commanded  rudder  6c'  independent  of  the 
disturbances  w. 

The  RMS  response  of  the  Tokyo  Maru  at  Fn=0.103  and  H/T=1.89  to  the  modeled 
design  disturbances  using  measurements  of  ip',  r',  and  n'  and  the  above 
control  and  filter  gains  is  shown  in  the  middle  column  of  Table  11.  The  RMS 
cost  established  using  eq.  (72)  is  also  shown.  Dimensionally  this  response 
is  $=0.43°,  ij=3.55m,  and  iS=2.6°.  The  optimal  stochastic  control  system 
is  therefore  very  effective  in  controlling  the  ship  when  subjected  to  the 
modeled  design  disturbances.  Notice  that  by  properly  processing  the  three 
available  measurements  the  controller  can  control  n to  3.55m  even  though 
the  RMS  error  in  the  n measurement  is  10m.  Also  shown  in  Table  11  are  the 
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RMS  responses  of  the  Tokyo  Maru  to  the  modeled  design  disturbance  when  controlled 
by  optimal  stochastic  controllers  designed  using  the  same  A and  B matrices  but 
one  measurement  n'  or  four  measurements  (ip  ' , r ' , 8 1 , n ' ) . As  expected  from 
the  discussion  above  concerning  observability,  control  based  on  an  n' 
measurement  alone  is  very  ineffective;  i.e.,  n=20m.  Also,  the  addition  of  a 
8'  measurement  produces  very  little  improvement  in  controller  effectiveness. 


measurements 

n'  only 

ip'  , r ' , n' 

ip',  r',  8' 

n' 

r 

0.04643 

0.00755 

0.00756 

r’ 

0.05145 

0.00876 

0.00876 

8' 

0.01436 

0.00501 

0.00502 

n' 

0.06921 

0.01223 

0.01211 

6' 

0.11841 

0.04521 

0.04500 

V 

0.14275 

0.05180 

0.05154 

J 

0.06473 

0.00566 

0.00560 

Table  11.  RMS  Response  of  Tokyo  Maru  to  Design  Disturbance 

using  Various  Measurements 


A brief  look  at  the  eigenvalues  associated  with  the  optimal  stochastic 
controller  for  the  Tokyo  Maru  at  Fn=0.103  and  H/T=1.89  is  of  value  at  this 
point.  Table  12  shows  the  open-loop  eigenvalues  for  the  ship  in  this 
condition  as  presented  in  Section  2.  The  dominant  eigenvalue  14=0. 0237 
indicates  that  the  ship  is  course  unstable.  Table  12  also  shows  the  closed- 
loop  eigenvalues  for  the  ship  controlled  with  the  optimal  control  gains  C 
presented  above.  The  optimal  controller  design  guarantees  a stable  controlled 
ship  which  is  verified  by  these  results.  The  dominant  poles  have  a real 
part  of  -0.5538  which  indicates  good  course  stability.  For  the  Kalman-Bucy 
filter  to  be  effective  in  processing  the  noisy  measurements  to  produce  the 
state  estimate  x it  is  necessary  that  the  error  in  this  estimate  have 
dynamics  which  are  faster  than  the  ship.  If  this  is  true,  the  filter  will 
produce  an  acceptable  estimate  x quickly  enough  after  some  change  that  it 
can  be  used  effectively  to  control  the  ship.  The  estimate  error  eigenvalues 
are  shown  in  Table  12  and  the  dominant  eigenvalues  indicate  that  the  filter 
will  respond  about  twice  as  fast  as  the  ship.  This  should  prove  acceptable 
as  will  be  confirmed  by  computer  simulation  results  presented  below. 


open-loop 


closed-loop 


estimate  error 


0.0 

-.5538+  .8669  j 

-.9111+  1.157  j 

0.0 

0.0237 

-0.9534 

-1.824 

-2.317 

-2.311 

-2.638 

-5.281 

-7.470 

-5.281 

-1.000  shaping 

-1.000  shaping 

-7.208+  7.274  j 

-1.000  filters 

-1.000  filters 

Table  12.  Eigenvalues  for  Optimal  Design  for  Tokyo  Maru 

The  optimal  stochastic  path  controllers  as  developed  here  provide 
effective  path  control  provided  sufficient  measurements,  such  as  \p' , r', 
and  r)'»  are  utilized.  In  the  following  sections,  we  study  the  effect  of 
the  changes  in  ship  characteristics  due  to  changing  water  depth  and  ship 
speed  on  the  design  and  performance  of  these  controllers.  The  response  of 
these  controllers  to  initial  condition  errors  and  the  specific  disturbances 
shown  in  Fig.  4 as  obtained  by  computer  simulation  will  also  be  presented. 


. Effects  of  Speed  and  Water  Depth  on  Controller  Design  and  Performance 


1 


■! 

As  shown  in  Section  2,  the  characteristics  f and  y. . for  a surface 
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ship  change  considerably  as  vessel  speed  and  water  depth  change.  The  purpose 
of  this  section  is  to  establish  how  these  changes  affect  the  design  of  the 
optimal  stochastic  path  controllers  and  more  importantly  how  these  changes 
affect  the  performance  of  the  optimal  stochastic  path  controllers.  These 
results  help  establish  which  speed  and  depth  conditions  should  be  used  in 
design  if  constant  control  and  filter  gain  matrices  are  to  be  utilized. 

These  results  can  also  indicate  whether  or  not  an  adaptive  path  control 
system  is  necessary.  Three  measurements  ip' , r',  and  g ' are  used  through- 
out this  section. 

4 1 Tokyo  Maru  Controller  Design  and  Performance 

1 2 

The  shallow  water  data  presented  by  Fugino  ' for  the  Tokyo  Maru  is 
only  complete  for  the  single  Froude  number  Fn=0.103.  Therefore,  we  could 
only  evaluate  the  effect  of  water  depth  on  the  design  and  performance  of 
the  Tokyo  Maru  controllers.  The  versatility  of  the  OPTSYS  program  allows 
a very  meaningful  study  to  be  conducted.  As  presented  in  Section  3 for 
H/T=l . 89 , the  optimal  controller  gain  matrix  C and  the  optimal  Kalmar.-Bucy 
filter  gain  matrix  K can  be  established  for  any  water  depth  and  the  RMS 
response  of  the  controlled  ship  to  the  modeled  design  disturbances  can  be 
evaluated.  It  is  also  possible  to  establish  the  RMS  response  of  the 
controlled  ship  at  any  water  depth  using  a specified  set  of  gain  matrices. 

This  allows  an  evaluation  of  the  performance  of  a controller  designed  for  one 
water  depth  when  operating  at  another  water  depth. 

Optimal  stochastic  path  controllers  were  designed  for  the  Tokyo  Maru 

for  the  five  depth-to-draft  ratios  for  which  specific  data  were  presented 

in  Section  2.  Table  13  gives  the  optimal  control  gains  as  a function  of 

H/T.  The  gain  on  n',  C.  , is  constant  with  respect  to  H/T  which  could 

simplify  implementation  of  an  adaptive  system.  Two  of  the  gains  on 

6'  and  C ^ on  Y'  show  major  changes  in  value  as  water  depth  changes. 

Gain  C,_  on  N ' shows  a somewhat  smaller  change.  Gains  C,  , , C._,  and  C,c 
Ib  11121b 

show  relatively  small  changes.  Table  14  gives  the  optimal  Kalman-Bucy  filter 
gains  for  the  three  water  depths  H/T=ao,  H/T=1.89,  and  H/T=1.30.  The 
gains  ^=0  for  all  depths  as  noted  in  Section  3 for  H/T=1.89.  In 
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general,  the  optimal  filter  gains  change  less  with  water  depth  than  do  the 
optimal  control  gains.  Some  gains  such  as  and  K^2  change  very  little; 
others  such  as  change  by  as  much  as  a factor  of  three.  The  most  important 
aspect,  however,  is  not  how  much  the  gains  should  change  with  water  depth  to 
remain  optimal  but  how  much  operating  at  an  incorrect  water  depth  affects 
the  control  performance. 


H/T 

1.30 

1.50 

1.89 

2.50 

CO 

C11 

5.605 

5.566 

5.647 

5.659 

5.984 

C12 

2.632 

2.722 

2.685 

2.561 

2.489 

C13 

10.985 

7.630 

6.315 

3.998 

1.887 

C14 

2.425 

2.425 

2.425 

2.425 

2.425 

C15 

-0.723 

-0.779 

-0.808 

-0.806 

-0.760 

C16 

547.1 

563.0 

675.1 

693.7 

731.6 

C17 

-79.25 

-56.05 

-50.15 

-26.29 

7.033 

Table  13.  Optimal  Control  Gains  for  Tokyo  Maru  versus  H/T 

The  important  considerations  of  how  water  depth  affects  controller 
performance  and  what  water  depth  should  be  used  in  a constant -gain  controller 
design  can  be  addressed  by  studying  Tables  15,  16,  and  17.  In  these  three 
Tables  we  present  the  RMS  response  of  optimal  stochastic  controllers  designed 
for  the  Tokyo  Maru  at  H/T=°°,  H/T=1.89,  and  H/T=1.30  when  the  ship  is 
subjected  to  the  modeled  design  disturbances  at  these  same  three  depths.  Table 
15  shows  the  RMS  response  and  RMS  cost  J when  the  ship  is  operating  at  H/T=® 
using  controllers  which  were  designed  to  be  optimal  at  H/T=1.30,  H/T=1.89,  and 
the  existing  H/T=°°.  Based  on  RMS  cost,  which  is  the  best  single  criterion, 
the  design  for  H/T=°°  is  superior  as  expected.  The  optimal  design  for  H/T=1.89 
shows  an  RMS  cost  increase  from  the  optimal  of  0.00681  to  0.00783  or  15  per 
cent.  The  RMS  lateral  offset  n increases  from  4.7  m to  6.5  m.  The  optimal 
design  for  H/T=l . 30  shows  a larger  RMS  cost  increase  to  0.01056  with  an  RMS 
lateral  offset  increase  to  8.95  m.  Table  16  shows  the  RMS  response  with  the 
same  three  control  systems  when  the  ship  is  operating  at  H/T=1.89.  The 
middle  column  for  the  optimal  design  is  a repeat  of  the  results  given  in 


optimal  filter  gains  K for  H/T=»J 


2.6319 

0.9158 

0.0013 

6.3892 

13.0320 

-0.0594 

1.5709 

2.9049 

-1.9634 

0.1245 

-0.8312 

1.9697 

0.0000 

0.0000 

0.0000 

0.0294 

0.1840 

0.0093 

-0.0180 

-0.1026 

0.0421 

optimal  filter  gains  K for  H/T=1.89 

2.6321 

0.9164 

0.0015 

6.3936 

12.2380 

-0.0724 

1.4384 

3.2622 

-1.2160 

0.1447 

-1.0125 

1.5367 

0.0000 

0.0000 

0.0000 

0.0263 

0.1673 

0.0130 

-0.0287 

-0.1710 

0.0319 

optimal  filter  gains  K for  H/T=1.30 

2.6252 

0.8906 

0.0037 

6.2137 

10.2795 

-0.0476 

0.8253 

1.9779 

-0.4600 

0.3619 

-0.6658 

0.9457 

0.0000 

0.0000 

0.0000 

0.0324 

0.1587 

0.0105 

-0.0406 

-0.1815 

0.0245 

Table  14.  Optimal  Filter  Gains  for  Tokyo  ktaru  versus  H/T 
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Section  3.  The  optimal  design  for  H/T=°°  shows  a decrease  in  the  RMS  lateral 
offset  n'  with  increased  iJj'  and  increased  rudder  usage.  Since  manipulation 
of  the  A and  B matrices  can  alter  the  relative  magnitudes  of  these  values,  this 
illustrates  the  need  to  utilize  a single,  rational  basis  upon  which  to  compare 
the  various  designs.  Based  on  the  RMS  cost,  the  design  for  H/T=°°  shows  an 
increase  from  the  optimal  0.00566  to  0.00644.  The  design  for  H/T=1.30  shows 
a smaller  increase  to  0.00611  with  an  increase  in  n'  but  decreased  ip'  and 
decreased  rudder  usage.  Table  17  shows  the  RMS  response  of  the  three  control 
systems  when  the  ship  is  operating  at  H/T=1.30.  The  design  for  H/T=1.89  shows 
a much  smaller  degredation  in  RMS  cost  at  this  shallower  depth  than  does  the 
design  for  H/T=°°.  Both  designs  show  improved  RMS  lateral  offset  g ' but 
increased  and  increased  rudder  usage. 


Design  for 
H/T=l . 30 

Design  for 
H/T=l . 89 

Optimal  Design 

r 

0.00936 

0.00850 

0.00965 

r' 

0.00749 

0.00753 

0.00982 

S' 

0.00718 

0.00756 

0.00834 

n' 

0.03087 

0.02238 

0.01607 

6 ' 

0.04548 

0.04533 

0.04773 

6 • 

c 

0.05288 

0.05249 

0.05406 

j 

0 . 01056 (+55%) 

0 . 00783 (+15%) 

0.00681 

Table  15 

RMS  Response  of 

Tokyo  Maru  Operating  at  H/T=°° 

Design  for 
H/T=l . 30 

Optimal  Design 

Design  for 
H/T=°° 

*' 

0.00675 

0.00755 

0.00989 

r ' 

0.00751 

0.00876 

0.01266 

6' 

0.00469 

0.00501 

0.00571 

n' 

0.01577 

0.01223 

0.01033 

6' 

0.04445 

0.04521 

0.05035 

V 

0.05124 

0.05180 

0.05639 

J 

0.00618 (+8%) 

0.00566 

0.00644 (+14*) 

Table  16.  RMS  Response  of  Tokyo  Maru  Operating  at  H/T=lii9 
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Optimal  Design 

Design  for 
H/T=l . 89 

Design  for 

H/T=°° 

r 

0.00697 

0.00790 

0.00904 

r' 

0.00830 

0.00970 

0.01228 

0.00234 

0.00248 

0.00274 

0.01133 

0.01083 

0.00984 

0.04318 

0.04459 

0.05017 

0.04850 

0.05003 

0.05562 

J 

0.00502 

0.00524 (+4%) 

0 . 00626 (+25%) 

Table  17.  RMS  Response  of  Tokyo  Maru  Operating  at  H/T=l . 30 


A study  of  the  results  presented  in  Tables  15,  16,  and  17  points  the  way 
to  a choice  of  design  water  depth  if  a constant-gain  stochastic  control  system 
is  to  be  utilized.  We  conclude  that  if  the  design  is  to  be  optimal  at  a 
single  water  depth,  it  should  be  designed  to  be  optimal  at  an  intermediate 
water  depth  where  the  ship  is  least  course  stable.  For  the  Tokyo  Maru,  the 
design  which  is  optimal  at  H/T=1.89  provides  the  best  overall  performance. 

All  three  controllers  are  stable  in  all  three  conditions  but  the  design  for 
H/T=1.89  is  developed  to  be  optimal  when  the  ship  is  most  difficult  to 

i control  and  this  seems  to  provide  the  best  performance  over  the  range  of 

water  depths.  For  ships  of  this  type,  however,  we  feel  that  the  results  of 
Tables  15,  16,  and  17  show  the  justification  for  the  development  of  an  adaptive 
system  which  could  automatically  adjust  the  C and  K gain  matrices  to  be 
optimal  at  any  operating  H/T. 

4 . 2 Mariner  Controller  Design  and  Performance 

The  Mariner  is  course  stable  in  all  water  depths  so  the  effect  of  water 
depth  on  the  controller  design  and  performance  could  be  expected  to  be  less 
than  shown  above  for  the  tanker  Tokyo  Maru.  Since  Fujino's  data  for  the  Mariner 
included  two  speeds,  we  conducted  a study  to  see  how  both  vessel  speed  and 
water  depth  affected  controller  design  and  performance  for  this  type  of  ship. 

This  was  begun  by  developing  optimal  designs  for  the  Mariner  at  the  four  combina- 
tions of  Froude  numbers  0.0905  (7  knots  full-scale)  and  0.155  (12  knots  full- 
scale)  and  water  depth-to-draft  ratios  H/T=°°  and  H/T=1.50.  Table  18  gives  the 
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resulting  optimal  control  gains  for  the  four  conditions.  The  results  are 
similar  to  those  shown  in  Table  13  for  the  Tokyo  Maru.  The  gain  on  g’,  C^, 
is  a constant  2.425  for  all  conditions  and  thus  appears  independent  of  ship 
type,  vessel  speed,  and  depth-to-draft  ratio.  Gains  on  0'  and  C on  r'  show 
a strong  dependence  on  both  speed  and  H/T.  Gain  shows  a larger  speed 
dependence  than  depth-to-draft  ratio  dependence;  gain  C ^ shows  the  opposite 
behavior . 


Table  18.  Optimal  Control  Gains  for  Mariner  at  Four  Conditions. 

Table  19  gives  the  optimal  Kalman-Bucy  filter  gains  for  the  Mariner  in 
the  four  water  depth  and  vessel  speed  conditions  studied.  As  with  the  Tokyo 
Maru  the  optimal  filter  gains  show  relatively  less  variation  with  depth-to- 
draft  ratio  than  do  the  optimal  control  gains.  There  appears  to  be  a greater 
vessel  speed  dependence  in  the  optimal  filter  gains  but  the  conclusion  varies 
from  term  to  term  as  can  be  seen  by  study  of  Table  19. 

To  study  the  effects  of  vessel  speed  and  water  depth  on  the  performance 
of  the  optimal  stochastic  controller  for  the  Mariner,  we  used  the  design  for 
Fn=0.155  and  H/T=°°  as  the  base  design.  This  was  the  least  course  stable 
condition  for  this  ship  as  shown  by  the  open-loop  eigenvalues  given  in  Table 
5.  The  RMS  response  of  the  ship  to  the  modeled  design  disturbances  while 
using  this  control  system  design  was  then  established  at  Fn=0.155  and  H/T= 
1.50,  at  Fn=0.0905  and  H/T=°°,  and  at  Fn=0.0905  and  H/T=1.50.  Table  20  gives 
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optimal  filter  gains  K for  Fn=0.155  H/T=°° 

1.3069 

0.9836 

0.0008 

1.6806 

23.3700 

-0.0118 

0.5647 

5.9375 

-1.0798 

0.2459 

-2.1920 

1.4612 

0.0000 

0.0000 

0.0000 

0.0029 

0.2727 

0.0036 

-0.0035 

-0.1745 

0.0133 

optimal  filter  gains  K for  Fn=0.0905 

H/T=°° 

2.2369 

0.9559 

0.0007 

4.7926 

18.5577 

-0.0218 

1.1006 

4.2767 

-1.4516 

0.2114 

-1.3801 

1.6954 

0.0000 

0.0000 

0.0000 

0.0138 

0.2061 

0.0047 

-0.0100 

-0.1269 

0.0199 

optimal  filter  gains  K for  Fn=0.155  H/T=1.50 

1 . 3066 

0.9780 

0.0016 

1.6711 

20.2666 

-0.0065 

0.2712 

4.1025 

-0.2179 

0.5045 

-1.2028 

0.6561 

0.0000 

0.0000 

0.0000 

0.0037 

0.2632 

0.0026 

-0.0049 

-0.1894 

0.0083 

optimal  filter  gains  K for  Fn=0.0905 

H/T=l . 50 

2.2357 

0.9451 

0.0010 

4.7386 

16.7166 

-0.0206 

0.7675 

3.8652 

-0.6359 

0.3313 

-1.3004 

1.1166 

0.0000 

0.0000 

0.0000 

0.0160 

0.1947 

0.0051 

-0.0155 

-0.1648 

0.0149 

Table  19.  Optimal  Filter  Gains  for  Mariner  at  Four  Conditions 

the  RMS  response  of  the  base  design  at  Fn=0.155  and  the  two  depths  H/T=°°  and 
H/T=1.50.  The  response  with  the  optimal  design  for  the  shallower  depth  is 
also  shown  for  comparison.  The  RMS  response  at  H/T=°°  is  ij/=92°,  n=4.1m.,  and 
6=3.64°  which  is  comparable  to  the  Tokyo  Maru  response  even  though  our  constant 
design  disturbance  appears  relatively  larger  to  the  Mariner.  At  the  shallower 
water  depth  the  base  design  shows  a decreased  RMS  lateral  offset  n‘  but 
increased  tji'  and  increased  rudder  usage  compared  to  the  optimal  as  was 


Operating 

H/T 

1.50 

1.50 

CO 

Design 

Fn 

Optimal 

0.155 

Optimal 

Condition 

H/T 

Design 

oo 

Design 

r 

0.01104 

0.01241 

0.01606 

r’ 

0.01531 

0.01920 

0.01835 

6’ 

0.00421 

0.00479 

0.01271 

n* 

0.01643 

0.01239 

0.02566 

6' 

0.05078 

0.05838 

0.06356 

6 

c 

0.06549 

0.07822 

0.08263 

j 

0.00858 

0.01058  (+23%) 

0.01500 

Table  20.  RMS  Response  of  Mariner  Operating  at  Fn=0.155 


Operating 

H/T 

1.50 

1.50 

OO 

OO 

Design 

Fn 

Optimal 

0.155 

Optimal 

0.155 

Condition 

H/T 

Design 

OO 

Design 

OO 

r 

0.00910 

0.00850 

0.01366 

0.01243 

r* 

0.01194 

0.01155 

0.01504 

0.01276 

6’ 

0.00557 

0.00564 

0.01184 

0.01144 

n' 

0.01422 

0.01032 

0.02142 

0.02245 

S' 

0.04636 

0.04963 

0.05717 

0.05551 

6 

c 

0.05340 

0.06255 

0.06667 

0.06837 

j 

0.00627 

0 . 00707 (+13%) 

0.01060 

0.01087 (+2.5%) 

Table  21.  RMS  Response  of  Mariner  Operating  at  Fn=0.0905 


shown  with  the  Tokyo  Maru.  The  RMS  cost  increases  23%  compared  to  the  optimal. 
Table  21  gives  the  RMS  response  of  the  base  design  when  the  ship  is  subjected 
to  the  modeled  design  disturbances  while  operating  at  Fn=0.0905  and  water 
depths  H/T=°°  and  H/T=1.50.  The  RMS  response  of  the  ship  with  the  optimal 
control  system  design  is  shown  in  both  conditions  for  comparison.  The  change 
in  speed  in  deep  water  shows  only  a small  2.5  per  cent  degredation  in  RMS 
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cost  compared  to  the  optimum.  The  RMS  lateral  offset  n'  decreases  but  ip ' 
and  rudder  usage  increase.  Notice  that  these  effects  of  slowing  down  are 
opposite  to  the  effects  of  entering  shallow  water.  This  apparently  explains 
why  the  degradation  in  RMS  cost  compared  to  the  optimum  is  less  at  slow 
speed  in  shallow  water,  Fn=0.0905  and  H/T=1.50,  than  at  the  control  system 
design  speed  Fn=0  .155  in  the  same  depth  water. 

In  general,  the  performance  of  a constant-gain  stochastic  control  system 
for  the  Mariner  shows  a smaller  loss  of  performance  as  vessel  speed  and  water 
depth  change  than  shown  for  the  Tokyo  Mar u as  water  depth  changes.  An 

I adaptive  system  might  still  be  justified  but  the  dimensional  response  is  good 

for  all  vessel  speeds  and  water  depths;  i.e.,  the  worst  case  gives  ri=4.1m., 
ijj=.920,  and  5=3.64°.  A constant-gain  design  should  be  a very  effective 
controller  for  a ship  such  as  the  Mariner  which  is  course  stable  under  all 
conditions . 


i 


5 . Simulation  of  Optimal  Stochastic  Path  Controllers 

This  section  presents  a general  discussion  of  the  digital  computer 
simulation  of  optimal  stochastic  path  controllers.  We  also  present  simulation 
results  for  the  controller  designs  developed  above  for  the  Tokyo  Maru.  These 
results  are  used  to  evaluate  the  validity  of  using  first-order  shaping  filters 
to  model  the  yawing  moment  and  sway  force  disturbances  in  the  controller 
design.  The  section  closes  with  a comparison  of  the  effectiveness  of  the 
optimal  controller  designed  for  the  Tokyo  Maru  at  H/T=°°  with  the  control 
which  might  be  provided  by  a human  helmsman. 

The  OPTSYS  optimal  stochastic  controller  design  program  can  evaluate 
the  RMS  response  of  the  controlled  system  to  the  design  stochastic  disturbances 
represented  by  the  power  spectral  density  Q while  using  measurements  which 
are  contaminated  by  noise  represented  by  the  power  spectral  density  R.  These 
RMS  results  are  very  useful  in  evaluating  the  designs  but  often  they  are 
difficult  for  designers  to  fully  interpret.  To  more  completely  show  system 
effectiveness  there  is  also  a need  to  simulate  the  response  of  the  controlled 
system  to  initial  condition  (start-up)  errors  and  specific  process  disturbances 
while  the  controller  is  receiving  measurements  containing  randomly  generated 
noise.  To  perform  the  simulations  needed  as  part  of  this  work,  we  developed 
the  SHIPSIM/OPTSIM  simulation  program  for  stationary,  linear,  optimal 
stochastic  control  systems. ^ This  program  uses  the  same  data  sets  used 
with  or  produced  by  the  OPTSYS  program  so  the  two  companion  programs  provide 
a valuable  capability  for  the  efficient  design  and  evaluation  of  these 
controllers. 

5 . 1 Formulation  of  Simulation  Equations 

The  physical  system  being  controlled  is  represented  by  the  system 
equations  (29);  i.e., 

x = FgX  + Gsu  + rsw  » x(tQ)  = xQ  , (73) 

and  the  measurement  equation  (36) , 

z = Hsx  + v . (74) 

The  optimal  stochastic  controller  design  is  actually  based  on  the  augmented 
or  est imator-design  state  equations  (42);  i.e.. 
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x'  = Hex’  + Geu  + Tew  , (75) 

whenever  shaping  filters  are  used  to  model  the  process  disturbances.  If 
shaping  filters  are  not  used  x'=x,  Fe=Fs,  Ge=Gs,  and  re=rs;  i.e.,  eq.  (73) 
and  eq.  (75)  are  identical.  The  optimal  stochastic  controller  design 
produces  the  optimal  control  law  eq.  (45);  i.e., 

u = Cx  , (76) 


and  the  Kalman-Bucy  filter  given  by  eq.  (50), 

k = Fex  + Geu  + K(z-Hex)  , x(t0)=x0  . (77) 

Substituting  eq.  (74)  and  eq.  (76)  into  eq.  (73)  and  eq.  (77)  produces  the 
state  and  state  estimate  equations  which  must  be  simulated. 


x = Fsx  + GsCx  + r sw  , x(t0)=Xo  , (78) 

x = Fe8  + GeCx  +KHsx  + Kv-KHex  , x(tQ)=x0  . (79) 

Equations  (78)  and  (79)  combine  to  give  the  (ns+ne)  system  of  coupled 
first-order  differential  equations  which  the  SHIPSIM/0PTS1M  program  simulates; 
i.e.  , 


(ns+ne) xl 


(q+p) xl 


d x 
dt'  x 


Fs  GsC 

r h 

X 

”rs 

o" 

w 

KHS  Fe+GeC-KHe 

X 

+ 

0 

K 

V 

(80) 


These  equations  are  integrated  by  SHIPS1M/OPTSIM  from  a user-specified 
initial  condition. 


X (to)"1 

II 

£o 

X (tQ) 

io 

(81) 


using  any  process  disturbance  vector  w(t)  produced  by  a user-supplied 
process  disturbance  subroutine.  The  zero-mean,  gaussian  measurement  noise 
vector  v(t)  is  randomly  generated  in  accordance  with  a user-specified 
standard  deviation  o_'  . SHIPSIM/OPTSIM  can  also  integrate  additional  differen- 
tial equations  and  calculate  any  non-integrated  quantities  desiied  as  part  of 
the  printed  or  plotted  output.  User's  Documentation,  Programmer's  Documentation, 
and  a listing  for  the  SHIPSIM/OPTSIM  program  have  been  published  separately.11 
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5 . 2 Selection  of  Integration  Method  and  Noise  Modeling 

For  a proper  simulation  of  eq.  (80) , the  selection  of  the  integration 
method  and  the  treatment  of  the  measurement  noise  must  be  performed  with 
care  to  ensure  that  the  correct  stochastic  response  is  achieved.  The  most 
direct  approach  is  to  utilize  a fixed  step-size  Euler  or  rectangular  inte- 
gration scheme.  This  has  the  effect  of  approximating  the  continuous 
Gauss-Markov  process  eq.  (80)  by  a discrete  Gauss-Markov  sequence.  In 
developing  the  design  of  the  optimal  control  system,  eq.  (66)  was  used  to 
establish  the  measurement  noise  power  spectral  densities.  The  noise  was 
first  assumed  to  be  exponentially  correlated  with  an  RMS  magnitude  Oj 
and  a correlation  time  Tj  which  is  much  faster  than  the  time  constants 
for  the  ship.  Equation  (66);  i.e., 

rjj  " 2 (° j 2 > x j , (82) 

was  then  used  to  calculate  the  magnitude  of  the  noise  power  spectral  density- 
This  was  then  assumed  constant  when  the  noise  was  approximated  as  white 
noise.  The  RMS  magnitude  Oj  used  in  eq.  (82)  is  not,  however,  the  correct 
standard  deviation  component  Oj'  to  use  in  the  simulation  of  the  system 
using  Euler  integration. 

The  use  of  a fixed  step-size  Euler  integration  has  the  effect  of 
approximating  the  continuous  Gauss-Markov  process  eq.  (73)  and  eq.  (74)  by 
the  discrete  Gauss-Markov  process  given  by. 


*i+l 

= *i2Li  + riwi  , 

i=0, 1 N 

(83) 

= H . x • + v.  , 

l—i  —l 

i=0,l, . . . ,N 

(84) 

E = Ci6ij  * (85) 

E^iVjT)  = Ri6ij  * (86) 

El”iZjTj  = E[!^i  <Xo-?o)T1  = E(vi  (?o"Fo)TI  = 0 » (87) 

where  Q,  and  R.  are  now  covariance  matrices  and  6. . is  the  Kronecker 
ii  i] 

delta  function.  The  control  can  be  omitted  in  this  discussion  without  loss 
of  generality.  If  we  take  as  the  error  covariance  of  the  discrete  Kalman 

filter  estimate  at  point  i,  the  covariance  can  be  propagated  to  the  next 
measurement  point  i+1  using  a time  update  to  give  the  error  covariance 
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prior  to  the  measurements  and  then  the  measurements  can  be  processed 

in  the  measurement  update  to  give  + ^ • Governing  update  equations  can 

be  expressed  as. 


T T 

W = Wi  + WA 

(time  update)  t 

(88) 

i+1  = Mi+1  - Ki+lHi+lMi+l 

(measurement  update)  , 

(89) 

where , 

Ki  = PiHiRi~1  = MiHj(HiMiH^+Ri)-1  . (90) 

When  Euler  integration  is  used  with  an  integration  step-size  At, 
x(t+At)  is  approximated  using  eq . (73)  to  be, 

x(ttAt)  = x(t)  + Fsx(t)At  + T w(t)At  , (91) 

— S— 

where  again  the  control  is  omitted  for  the  purposes  of  this  discussion.  The 
error  covariance  equation  (52)  is  not  included  in  the  simulation  but  the 
use  of  Euler  integration  has  the  effect  of  approximating  the  error  covariance 
P(t+At)  as, 

P (t+At)  = P(t)  + FP  (t ) At  + P (t)  FTAt  + rQl'TAt-K(t)HP(t)  At  , (92) 

where  the  filter  gain  matrix  is  given  by, 

K = P(t)HTR-1  . (93) 

For  a simulation  to  be  valid  it  is  necessary  that  the  discrete  error  covariance 
given  by  eq.  (88),  (89),  and  (90)  be  the  same  at  the  end  of  each  integration 
step  as  the  continuous  error  covariance  given  by  eq.  (92)  and  (93).  We  show 
in  ref.  11  that  this  is  true  if  the  disturbance  and  noise  covariance  matrices 
in  eq.  (85)  and  (86)  are  related  to  the  continuous  system  power  spectral 
density  matrices  as  used  in  the  system  design  by, 

QjAt  = Q , (94) 

and 

R. At  = R . (95) 

l 

The  correct  measurement  noise  covariance  to  use  in  a simulation  thus 

depends  on  the  integration  step-size  At.  If  o_.’  is  the  standard  deviation 
for  noise  element  j (square  root  of  diagonal  element  j of  R^  which  is 


with  the  integration  performed  by  Euler  integration  using  a fixed  step-size 
At.  SHIPSIM/OPTSIM  uses  a fixed  step-size  Euler  integration  and  in  our 
simulations,  eq,  (96)  was  used  to  relate  the  noise  standard  deviation  vector 
a_'  to  the  design  measurement  noise  level  represented  by  the  power  spectral 
density  R.  A non-dimensional  step-size  of  At=0.005  was  found  to  be 
satisfactory  for  most  of  these  simulations. 


The  SHIPSIM/OPTSIM  program  was  used  to  perform  a number  of  simulations 
of  the  Tokyo  Mar u controlled  by  the  optimal  controllers  we  have  developed. 

Some  of  the  results  for  the  Tokyo  Maru  operating  at  a water  depth-to-draft 
ratio  H/T=1.89  are  presented  here  to  further  illustrate  the  effectiveness 
of  these  controllers.  Recall  that  the  controllers  we  have  studied  are  for 
stationary  systems;  i.e.,  the  statistics  of  the  disturbance  and  measurement 
noise  (RandQ)  and  the  various  system  matrices  are  assumed  constant.  The 
resulting  control  feedback  gains  and  Kalman-Bucy  filter  gains  are  also 
constant . 

An  important  question  with  a stationary,  constant-gain  control  system 

is  how  the  controlled  system  will  respond  if  the  controller  is  first  turned 

on  while  an  initial  condition  error  exists.  We  simulated  the  start-up  response 

of  the  Tokyo  Maru  controlled  by  the  optimal  stochastic  controller  designed 

for  the  existing  depth  H/T=1.89  as  presented  in  Section  3.  All  initial  states 

were  taken  as  zero  except  the  lateral  offset  n which  was  taken  as  one  half 

o 

breadth  B/2  (23.75  m.  full-scale).  The  filter  has  no  prior  knowledge  at 
start-up  so  the  estimate  of  all  the  states  x was  taken  as  zero.  This 
corresponds  to  the  situation  where  the  ship  is  just  entering  the  restricted 
path  area  and  is  on  a path  parallel  to  but  offset  B/2  to  starboard  (r^O)  of 
the  prescribed  path  when  the  controller  is  turned  on.  The  lateral  offset 
response  and  rudder  angle  response  are  shown  in  Fig.  6 and  Fig.  7,  resoectively . 
The  lateral  offset  is  unchanged  for  about  one  ship  lenqth  while  the  filter 
generates  an  estimate  of  the  initial  lateral  offset  error,  the  rudder  responds 


to  a commanded  turn  to  port,  and  then  the  ship  responds.  The  ship  crosses 
the  prescribed  path  after  four  ship  lengths  (t'=4.0)  and  then  overshoots 
only  about  3.2  m.  before  settling  onto  the  prescribed  path.  The  small 
oscillations  about  the  path  after  about  t'=10  are  due  primarily  to  the  measure 
ment  noise  which  is  present  throughout  the  simulation.  The  excellent  startup 
response  with  the  optimal  controller  is  clearly  shown  by  Fig.  6.  The  maximum 
rudder  usage  shown  on  Fig.  7 is  dimensionally  about  3.4°.  The  continuing 
rudder  movement  due  to  the  measurement  noise  after  about  t'=10  has  peak  values 
of  only  about  1°  and  while  this  would  be  undesirable  for  steering  system  wear 
reasons  in  continuous  operation  it  should  be  very  acceptable  for  short-term 
operations  in  restricted  waters. 


T -N.  NOISE.  89.  C4K  FOR  H/T-1.89.  INITIAL  OFFSET  O.S  BEAN 


Figure  6.  Lateral  Offset  Response  to  B/2  Initial  Offset:  Optimal  Control 


The  design  process  disturbance  used  to  develop  the  optimal  controllers 
was  based  on  the  yaw  moment  and  sway  force  due  to  a passing  ship  as  shown  in 
Fig.  4.  These  disturbances  were  modeled  by  two  first-order  shaping  filters 
as  described  in  Section  3.  Since  it  is  of  interest  to  see  how  the  optimally 
controlled  system  would  respond  to  the  actual  yaw  moment  and  sway  force 
disturbances  shown  in  Fig.  4,  the  Tokyo  Maru  with  the  optimal  controller 
was  next  simulated  to  pass  the  ship  respresented  by  these  disturbances.  The 
simulation  was  begun  with  no  initial  error;  i.e.,  x=ic=0,  and  the  two  ships 
were  simulated  to  passed  beam-to-beam  at  t'  = 7.  Simulation  results  are  shown 
in  Fig.  8,  9,  and  10.  Figure  8 shows  the  lateral  offset  response  of  the 
Tokyo  Maru.  The  maximum  response  is  only  about  1.7m.  which  occurs  at  about 
t'=8.5  or  one  ship  length  after  the  ship  has  passed.  This  response  is  low 
for  two  reasons.  First,  the  disturbances  shown  in  Fig.  4 are  relatively 
small  for  a ship  the  size  of  the  Tokyo  Maru.  Secondly,  the  optimal  control 
system  is  very  effective  in  controlling  the  ship.  While  response  to  the 
measurement  noise  after  t'-18  appears  large  due  to  the  plot  scaling,  it  has 
a peak  value  of  only  about  0.6  m.  The  maximum  response  due  to  the  passing 
ship  would  increase  roughly  linearly  with  the  disturbance  magnitude  if  the 
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passing  ship  were  larger  or  closer  or  if  the  effect  of  shallow  water  were  to 

increase  the  disturbance  magnitude  significantly  as  discussed  in  Section 

3.2.  The  disturbance  magnitudes  shown  in  Fig.  4,  are,  therefore,  completely 

valid  for  linear  control  system  design  and  evaluation.  The  rudder  angle 

response  is  shown  in  Fig.  9.  The  maximum  rudder  angle  of  about  6.9°  occurs 

just  a few  seconds  after  the  maximum  disturbance  at t'  =7  when  the  ships  are 

passing  beam-to-beam.  The  estimate  of  the  yaw  moment  disturbance  N produced 

by  the  Kalman-Bucy  filter  is  shown  in  Fig.  10.  Notice  that  the  maximum 
-5 

value  of  -19.4x10  at  t'=7  compares  very  favorably  with  the  peak  value  of 
-20.0x10  ^ shown  on  Fig.  4.  This  peak  value  also  occurs  when  the  ships  are 
beam-to-beam  or  t'=7  in  this  simulation.  The  estimate  of  the  sway  force 
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T-H.N0IX.H/T-1.M.C4*  FOB  h/T- 1.88.  PASSING  SHIP  DISTURBANCE 

Figure  9.  Rudder  Angle  Response  to  Passing  Ship:  Optimal  Control 


TIME 

T-fl.N0IX.NFT«t.M.C4K  FOR  H/1-1. 88. PASSING  SHIP  01STIMANCE 


Figure  10.  Yaw  Moment  Estimate  for  Passing  Ship:  Optimal  Control 
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The  particular  lateral  offset  response  shown  in  Fig.  8 is  strongly 
affected  by  the  sequence  of  randomly  generated  measurement  noise  included  in 
our  simulation.  SHIPSIM/OPTSIM^  uses  the  IBM  Scientific  Subrouting  Package 
random  number  generator  subroutine  RANDU  to  generate  the  sequence  of  random 
numbers  which  establish  the  noise  level  for  each  measurement  at  each  integration 
step  in  accordance  with  the  user-specified  standard  deviation  o_'  . To  allow 
a valid  comparison  among  competing  control  systems  SHIPSIM/OPTSIM  always  begins 
with  the  same  initial  random  number  seed.  Thus,  all  runs  utilize  the  same 
sequence  of  randomly  generated  numbers  unless  the  program  is  specifically 
modified  to  change  the  initial  seed.  When  the  disturbance  magnitude  is 
fairly  small  as  used  here,  the  response  due  to  the  disturbance  is  roughly 
comparable  to  the  response  due  to  the  measurement  noise  and  the  maximum 
system  response  is  noticeably  affected  by  the  particular  sequence  of  measurement 
noise  used  in  the  simulation.  The  response  shown  in  Fig.  8 is  such  that  the 
ship  is  already  moving  to  port  (g<0)  due  to  measurement  noise  when  the 
disturbnace  due  to  the  passing  ship  is  first  felt  (t'=5).  This  causes  the 
maximum  lateral  offset  of  the  ship,  which  is  also  to  port,  to  be  increased. 

To  clarify  this.  Fig.  11  shows  the  lateral  offset  response  of  the  Tokyo  Maru 
to  the  same  passing  ship  disturbance  when  there  is  no  noise  in  the  measure- 
ments received  by  the  filter;  i.e.,  a_'=0.  The  maximum  lateral  offset  is  about 
1.5  m.  in  lieu  of  1.7  m.  as  shown  in  Fig.  8.  It  is  also  possible  for  the 
measurement  noise  to  reduce  the  maximum  lateral  offset.  This  is  illustrated 
by  the  lateral  offset  response  in  the  simulation  shown  in  Fig.  12.  This 
simulation  was  conducted  with  a modified  initial  random  number  seed.  For 
this  sequence  of  measurement  noise,  the  ship  happens  to  be  moving  to  star- 
board (n>0)  due  to  the  measurement  noise  when  the  disturbance  due  to  the 
passing  ship  is  first  felt.  This  causes  the  maximum  lateral  offset  to  port 
due  to  the  passing  ship  to  be  reduced  to  only  about  1.2  m.  compared  to  the 
1.5  m.  shown  in  Fig.  II  with  no  measurement  noise  and  the  1.7  m.  shown  in 
Fig.  8. 

Simulation  results  can  also  show  the  degredation  in  controlled  system 
response  which  would  occur  if  the  controller  were  designed  to  be  optimal  at 
some  water  depth  other  than  the  water  depth  in  which  the  ship  is  operating. 

We  simulated  the  response  of  the  Tokyo  Maru  operating  at  a water  depth-to-draf t 
ratio  H/T=1.89  while  being  controlled  by  a controller  which  was  designed  to 
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to  be  optimal  in  deep  water  H/T=°°.  Figure  13  shows  the  lateral  offset 
response  to  a B/2  initial  lateral  offset  error.  This  response  can  be 
compared  directly  with  the  response  shown  in  Fig.  6 for  the  Tokyo  Maru  when 
the  optimal  control  and  filter  gains  are  used.  The  initial  return  to  the 
prescribed  path  is  noticeably  faster  with  the  non-optimal,  deep  water 
controller.  The  maximum  overshoot  however,  is,  about  15.4  m.  in  lieu  of  the 
overshoot  of  only  about  3.2  m.  shown  in  Fig.  6 and  it  takes  the  ship  about 
4 ship  lengths  longer  to  settle  onto  the  prescribed  path.  The  loss  in 
performance  with  the  use  of  non-optimal  control  and  filter  gains  at  H/T=1.89 
is  evident  from  a comparison  of  these  two  figures.  Figure  14  shows  the  lateral 
offset  response  to  a passing  ship  as  represented  by  the  yaw  moment  and  sway 
force  disturbance  shown  in  Fig.  4 when  the  Tokyo  Maru  is  again  controlled  by 
a controller  which  was  designed  to  be  optimal  at  H/T=°°.  This  response  can  be 
compared  directly  with  the  response  shown  in  Fig.  8 for  the  ToKyo  Maru  when 
controlled  by  the  optimal  controller  for  the  existing  H/T=1.89  water  depth. 

In  this  situation,  the  maximum  response  does  not  degrade  as  much  as  in  the 
initial  offset  simulation.  The  maximum  lateral  offset  is  about  1.9  m.  with 
the  non-optimal  H/T=°°  design  compared  to  about  1.7  m.  with  the  optimal  H/T=1.89 
design  controller.  This  result  could  be  expected  based  on  the  RMS  response 
results  for  the  Tokyo  Maru  controlled  by  these  two  controllers  at  H/T=1.89  as 
shown  in  Table  16.  There  the  design  for  H/T=“>  showed  a 14  percent  increase 
in  RMS  cost  but  a 16  percent  decrease  in  RMS  lateral  offset  compared  to  the 
controller  designed  to  be  optimal  at  H/T=1.89. 

In  general,  we  feel  that  these  simulation  results  confirm  that  the 
optimal  stochastic  path  controllers  developed  here  would  provide  very 
effective  ship  path  control  in  restricted  waters.  The  effectiveness  of  these 
controllers  will  be  further  compared  to  the  control  which  could  be  expected 
from  a human  helmsman  in  Section  5.5  below. 
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Figure  13.  Lateral  Offset  Response  to  B/2  Initial  Offset:  h/T=°°  Optimal  Design 
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Figure  14.  Lateral  Offset  Response  to  Passing  Ship:  H/T~™  Optimal  Design 
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5 . 4 Verification  of  Disturbance  Modeling 

In  this  section,  we  look  a little  more  closely  at  the  effectiveness  of  the 
Kalman-Bucy  filter  to  estimate  the  state  of  the  ship.  In  particular,  the 
effectiveness  of  the  filter  to  estimate  the  yaw  moment  and  sway  force  acting 
on  the  ship  due  to  a passing  ship  is  studied  in  more  detail  as  a means  of 
verifying  the  validity  of  using  first-order  shaping  filters  to  model  the  yaw 
moment  and  sway  force  disturbances.  We  also  compare  the  effectiveness  of 
the  optimal  controller  design  presented  in  Section  3 with  that  of  an  optimal 
controller  designed  without  the  use  of  shaping  filters;  i.e.,  the  yaw  moment 
and  sway  force  disturbances  are  assumed  to  be  white  noise  instead  of  exponential 
correlated  disturbances.  We  feel  these  results  clearly  show  the  validity 
of  our  disturbance  modeling  approach. 

Figure  15  shows  a more  detailed  presentation  of  the  yaw  moment  distur- 
bance estimates  N'  produced  by  the  optimal  Kalman-Bucy  filter  in  the 
simulation  runs  shown  in  Fig.  8,  Fig.  11,  and  Fig.  12.  Also  shown  in  Fig. 

15  for  comparison  is  the  yaw  moment  disturbance  N'  actually  used  in  these 
simulations  as  shown  in  Fig.  4.  The  yaw  moment  disturbance  estimate  for  the 
simulation  run  without  measurement  noise;  i.e.,  o'=0,  represents  the  expected 

average  of  estimates  over  a large  number  of  runs.  This  estimate  was  obtained 
in  the  simulation  run  shown  in  Fig.  11.  Comparing  this  estimate  with  the 
actual  yaw  moment  disturbance  shows  the  filter  to  be  generally  very  effective 
in  estimating  the  yaw  moment  disturbance.  This  is  remarkably  fast  when  one 
considers  that  the  filter  can  only  estimate  N’  after  the  ship  begins  to 
respond  to  the  disturbance  and  this  is  reflected  in  the  measurements  ijj ' , r', 
and  n'-  The  estimate  is  low  in  peak  value  by  about  20%  for  the  initial 
positive  peak  and  by  about  13%  for  the  maximum  disturbance  when  the  ships  are 
beam-to-beam  at  t'=7.  The  filter  is  less  effective  at  the  final  positive 
peak  where  the  estimate  is  low  by  about  a factor  of  4.  The  yaw  moment 
disturbance  estimates  for  the  simulation  run  using  the  standard  measurement 
noise  initial  random  number  seed  as  shown  in  Fig.  8 and  for  the  simulation 
run  using  the  modified  seed  as  shown  in  Fig.  12  are  also  shown  in  Fig.  15  to 
illustrate  the  effect  of  specific  sequences  of  measurement  noise  on  the 
estimate  N'  produced  by  the  Kalman-Bucy  filter. 
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Figure  16  shows  the  sway  force  disturbance  estimates  Y'  produced  by 
the  optimal  Kalman-Bucy  filter  in  the  simulation  runs  shown  in  Fig.  8, 

Fig.  11,  and  Fig.  12.  The  actual  sway  force  disturbance  Y'  used  in  these 
simulations  as  shown  in  Fig.  4 is  also  shown  for  comparison.  Again  the 
sway  force  disturbance  estimate  for  the  simulation  run  without  measurement 
noise  represents  the  expected  average  of  estimates  over  a large  number  of 
runs.  In  this  situation,  the  filter  is  much  less  effective  in  estimating  the 
sway  force  disturbance  than  it  is  in  estimating  the  yaw  moment  disturbance. 

The  filter  is  about  0.2  ship  lengths  behind  and  about  a factor  of  3 low  in 
estimating  the  peak  sway  force  disturbance  which  occurs  when  the  ships  are 
beam-to-beam.  Since  the  optimal  stochastic  controller  provides  excellent 
response  in  general,  this  inability  to  estimate  Y'  as  accurately  as  N'  does 
not  seem  to  be  a practical  problem. 

The  lack  of  ability  of  the  filter  to  effectively  estimate  Y'  is  at 
first  not  apparent  from  a look  at  the  observability  of  each  mode  rf  the  system 
with  each  measurement  as  shown  in  Table  9.  This  is  because  each  row  in 
Table  9 reflects  the  relative  observability  of  each  mode  with  each  measurement. 
If  this  Table  included  measurements  of  all  seven  states,  the  elements  of  each 


50 


-56- 


Figure  16.  Sway  Force  Disturbance  Estimates  for  Passing  Ship 

n 

row  would  be  seen  to  be  orthonormal ized;  i.e.  . E cos^  8. . =1  , when  the 

1=0  il 

measurements  are  just  the  states.  These  data  do  not  reflect  the  absolute 
observability  of  one  mode  compared  to  another.  Comparisons  made  within  each 
row  of  Table  9 are  therefore  valid  but  numerical  comparisons  within  each 
column  of  Table  9 are  not  valid.  The  primary  reason  for  the  poor  estimate 
of  Y'  is  due  to  speed  at  which  the  disturbance  occurs  compared  to  the 
time  constants  of  the  controlled  ship.  The  Y'  disturbance  occurs  about 
twice  as  fast  as  the  N'  disturbance  and  as  a result  is  too  fast  to  have 
much  of  an  effect  on  the  ship.  This  Y'  disturbance  therefore  has  a low 
level  of  absolute  observability.  A test  run  using  a Y'  disturbance  which 
is  identical  to  the  N'  disturbance  shown  in  Fig.  4 resulted  in  a good 
estimate  Y'  and  confirms  this  conclusion.  For  the  same  physical  reasons 
that  Y'  in  Fig.  4 does  not  affect  r',  ip ' , and  n * significantly  resulting 
in  a low  level  of  absolute  observability,  a weak  estimate  Y'  does  not 
seriously  affect  the  overall  controller  performance. 

The  effectiveness  and  validity  of  modeling  yaw  moment  and  sway  force 


disturbances  as  first-order  filters  in  the  optimal  stochastic  controller 
design  can  be  further  illustrated  by  comparing  this  controller  with  an  optimal 
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controller  design  developed  for  the  Tokyo  Maru  at  H/T=1.89  without  the  use 
of  the  shaping  filters.  In  this  case,  the  yaw  moment  and  sway  force  distur- 
bances were  considered  to  be  white  noise  with  power  spectral  densities  given 
by  eq.  (56)  and  eq.  (57) . The  resulting  optimal  design  produces  feedback 
gains  on  ip' , r‘,  6',  n',  and  6'  which  are  identical  to  those  presented  in 
Section  3 for  the  design  which  utilizes  the  shaping  filter  modeling  of  the 
disturbances.  The  closed-loop  eigenvalues  are  therefore  identical  to  the 
first  five  eigenvalues  shown  for  this  design  in  Table  12.  The  diim  ision  of 
the  Kalman-Bucy  filter  is  reduced  to  5;  i.e.,  x includes  just  the  system 
states  and  not  the  yaw  moment  and  sway  force  disturbances.  Table  22  shows 
the  RMS  response  of  the  optimal  controller  designed  for  the  Tokyo  Maru  at 
H/T=1.89  without  shaping  filters  with  that  of  the  optimal  controller  designed 
using  shaping  filters.  This  latter  data  is  a repeat  of  that  given  in  Tables 
11  and  16.  The  RMS  cost  increases  3.5  times  when  the  shaping  filters  are  not 
used  in  the  design;  the  RMS  lateral  offset  from  the  prescribed  path  increases 
from  3.55  m.  when  shaping  filters  are  used  to  8.61  m.  The  control  system 
developed  without  the  use  of  the  shaping  filters  is  much  less  effective 
based  on  the  RMS  response  to  the  design  disturbances. 


without  shaping  filters 

with  shaping  filters 

r 

.02108  (1.21°) 

.00755  (.43°) 

r' 

.03342 

.00876 

6' 

.00898 

.00501 

V 

.02968  (8.61m) 

.01223  (3.55m) 

6' 

.07651  (4.4°) 

.04521  (2.6°) 

<5  * 
c 

.09044 

.05180 

j 

.01966 

.00566 

Table  22.  RMS  Response  with  Optimal  Controllers  Designed 
with  and  without  Shaping  Filters 


A comparison  of  the  effectiveness  of  the  optimal  controllers  designed 
with  and  without  the  use  of  shaping  filters  to  model  the  yaw  moment  and  sway 
force  disturbances  is  perhaps  clearest  when  simulation  results  are  studied. 
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Figure  17  shows  the  lateral  offset  response  of  the  Tokyo  Maru  at  H/T=1.89 
to  the  passing  ship  disturbance  shown  in  Fig.  4 when  the  ship  is  controlled 
by  the  optimal  controller  developed  without  the  use  of  shaping  filters  to 
model  the  disturbances.  This  response  can  be  compared  directly  with  the 
response  shown  in  Fig.  8 for  the  Tokyo  Maru  controlled  by  the  optimal  controller 
developed  with  the  use  of  shaping  filters.  The  maximum  lateral  offset  in 
Fiq.  17  is  4.9  m.  compared  to  1.7  m.  shown  in  Fig.  8.  Thus  the  use  of 
shaping  filters  to  model  the  disturbances  results  in  an  optimal  stochastic 
controller  design  which  is  almost  three  times  as  effective  in  limiting  the 
lateral  offset  from  the  prescribed  path  when  a ship  passes.  Figure  18  shows 
the  rudder  usage  in  the  simulation  shown  in  Fig.  17.  The  maximum  rudder 
angle  is  about  5.3°  compared  to  a maximum  of  about  6.9°  shown  in  Fig.  9. 


Lateral  Offset  Response  to  Passing  Ship:  Optimal  Design 
without  Shap’ng  Filters 


Figure  17. 
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Figure  18.  Rudder  Angle  Response  to  Passing  Ship:  Optimal  Design 

without  Shaping  Filters 


5 . 5 Comparison  with  Human  Helmsman  Control 

The  final  simulation  results  to  be  presented  area  comparison  between  the 
control  provided  by  the  optimal  stochastic  controllers  developed  here  and  the 
control  which  might  be  expected  by  a human  helmsman.  To  be  judged  effective, 
the  stochastic  controllers  should  probably  be  at  least  as  good  as  a human. 

If  the  two  were  about  equal  based  on  typical  performance,  the  automated 
control  might  still  be  justified  and  desirable  since  it  is  the  "human  error" 
and  occasionally  bad,  non-typical  human  performance  which  usually  gets  a 
ship  into  trouble.  An  automated  system  could  provide  a more  consistent, 
repeatable  performance  which  would  be  less  affected  by  fatigue,  stress, 
distractions,  etc. 

In  order  to  perform  a digital  computer  simulation  representing  human 
helmsman  control,  a human  model  must  be  available.  Suitable'  human  models  for 
ship  control  research  ate  very  limited.  Human  models  for  the  ship  hradirui 

2 ri  2 (i 

control  problem  have  been  presented  by  Moot t and  Stuurman.  These  can  be 
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shown  to  be  essentially  equivalent.  Stuurman  presents  data  for  a single 
ship.  Hooft  presents  specific  data  for  helmsmen  control  for  a wide  range  of 
course-stable  ships.  For  our  work,  we  adapted  Hooft' s model  to  apply  to  path 
control  and  have  extrapolated  his  data  to  apply  to  the  Tokyo  Maru  in  deep 
water  where  the  ship  is  course  stable.  It  would  be  expected  that  the  helmsman 
would  be  less  effective  in  controlling  the  Tokyo  Mar u at  H/T=1.89  where  it  is 
course  unstable  but  the  human  model  data  necessary  to  allow  simulation  at 
this  water  depth  is  not  available. 

Hooft's  model  for  human  heading  control  was  derived  from  the  response  of 
helmsmen  who  were  instructed  to  steer  a ship  in  a land-based  simulator  along 
a prescribed  sinusoidal  path.  In  the  modeling,  he  first  assumed  that  the 
rudder  angle  6 would  be  positioned  proportional  to  the  course  error  A\J; 
and  the  yaw  rate  error  Ar;  i.e., 

6 = k1A0  + k2Ar  . (97) 

He  then  assumed  the  helmsman  would  compare  the  existing  rudder  angle  6 with 
that  required  rudder  angle  6 *=Jc needed  to  steer  along  the  prescribed 
sinusoidal  course.  This  gives  an  error  Detween  the  existing  and  required 
rudder  angles, 

AS*  = 5*-5  , (98) 


or, 


k -jAtJ;  = k3<J/-S 

which  gives. 

Alp  = iJ;-5/k3 


(99) 


(100) 


Combining  eq.  (97)  and  (100)  yields  Hooft's  helmsman  heading  control  transfer 
function. 


G(s) 


6(s) 

AiJMs) 


klk3 

kl+k3 


s+1 


kl+k3 


S+1 


Tls+1 
11  tNs+1 


(101) 


Hooft  conducted  the  simulator  experiments  with  helmsmen  to  establish  the 
parameters  a,  Tp»  and  TN.  Since  the  helmsmen  respond  differently 
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j 


when  controlling  different  types  of  ships,  Hooft  obtained  his  data  as  functions 

27 

of  the  K and  T in  Nomoto's  maneuvering  equation  ; i.e., 

Tr  + r = K6  . (102) 

The  data  a,  TL,  and  TN  therefore  apply  to  a wide  range  of  ship  dynamics  but 
were  unfortunately  limited  to  course-stable  ships. 

Hooft's  model  does  not  apply  directly  to  ship  path  control  but  provides 
the  best  available  helmsman  model.  We  therefore  adapted  the  model  to  path 
control  by  assuming  that  the  transfer  function  eq.  (101)  applied  not  just  to 
a heading  error  Ai p but  to  an  error  function  which  includes  a weighted  sum 
of  a heading  error,  a yaw  rate  error,  and  a lateral  deviation  error.  We 
thus  used, 

(103) 


G (s)  = 


6 (s) 
e (s) 


= a 


tls+1 

tns+i 


where  the  error  function  was  given  by, 

e = (tp-ip0)+E2  (r-r0)+B3(n-n0)  = B^+b^+B^  , (104) 

since  ^0=r0=r|0=®  on  prescribed  path.  To  implement  the  transfer  function 

in  the  simulation,  eq.  (103)  was  first  written. 


T s+’/T 

G (s ) = a — — r - = A 


s+1/Tt 


tn  s+Vtn  s+VTn 


(105) 


which  can  then  be  implemented  by  an  additional  state-variable  equation 
representing  an  internal  helmsman  state, 

y'  = -jrr  y*  + A(B1i(/,+B2r'+B3n’)  , 


(106) 


and  a rudder  command  equation, 


V = 


ml  — rP  * 

N L 


W 


y'  + A(B1i|/’+B2r '+B3n' ) 


(107) 


These  equations  were  used  in  a simulation  to  approximate  helmsman  path  control 
of  the  Tokyo  Maru  at  H/T=ao.  In  nondimensional  form,  Hooft's  data  yields  A=13.20, 
T^j=0.6630,  and  T^=0.5835  for  the  Tokyo  Maru  characteristics  at  H/T=“>. 
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Since  the  helmsman  model  was  only  available  for  a course-stable  condition, 
the  helmsman  results  cannot  be  compared  directly  with  the  optimal  stochastic 
controller  simulations  performed  at  H/T=1.89  and  presented  above.  The  optimal 
stochastic  controller  designed  for  deep  water  was  therefore  simulated  with 
the  passing  ship  disturbance  shown  in  Fig.  4 to  provide  a basis  for  comparison 
with  helmsman  control.  Lateral  offset  response  is  shown  in  Fig.  19.  Since 
the  Tokyo  Maru  is  course  stable  at  H/T=°°,  the  optimally  controlled  ship  is 
also  more  course  stable  than  at  H/T=1.89.  The  maximum  lateral  offset  is  only 
0.75  m.  compared  to  the  1.7  m.  shown  in  Fig.  8 for  H/T=1.89.  The  rudder 
angle  response  is  shown  in  Fig.  20.  The  maximum  rudder  angle  is  about  6.4°. 
The  degredation  of  controller  performance  with  changing  water  depth  can  be 
further  seen  by  comparing  the  results  shown  in  Fig.  19  and  Fig.  14.  The 
performance  of  the  optimal  controller  operating  at  its  design  depth  H/T=°° 
degrades  in  maximum  lateral  offset  response  from  the  0.75  m.  shown  in  Fig. 

19  to  the  maximum  response  of  1.9  m.  shown  in  Fig.  14  as  the  ship  is  brought 
into  the  shallower,  H/T=1.89  water. 


Figure  19.  Lateral  Offset  Response  to  Passing  Ship:  Optimal  Control  at  H/T=rn 
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The  response  of  the  Tokyo  Maru  to  the  passing  ship  disturbance  shown  in 
Fig.  4 while  operating  at  H/T=“>  under  control  of  the  helmsman  defined  by  eq. 

(106)  and  (107)  is  shown  in  Fig.  21  and  Fig.  22.  For  these  simulations  we 
have  used  B^=0.8,  B2=0.1,  and  B(=0.1  . These  quantities  include  both 
weighting  considerations  and  scaling  differences  among  ip‘,  r',  and  n 1 errors. 
For  these  simulations,  the  helmsmen  was  given  the  same  noisy  measurements 
used  by  the  stochastic  controller.  This  is  very  conservative  since  human 
reaction  to  a measurement  display  will  introduce  additional  error.  Further, 
the  helmsman  will  not  read  each  measurement  at  each  decision  point.  This 
will  be  particular  evident  in  the  n measurement  if  there  is  no  automated 
output  of  the  lateral  offset  from  the  prescribed  path  and  the  helmsman  or  an 
assistant  must  establish  n manually  from  a radar  screen.  The  lateral  offset 
response  in  Fig.  21  can  be  compared  directly  with  the  response  with  the 
optimal  stochastic  controller  as  shown  in  Fig.  19.  The  maximum  lateral  offset 
response  is  1.5  m.  compared  to  the  0.75  m.  shown  in  Fig.  19.  The  maximum 
rudder  angle  response  in  Fig.  22  is  about  10.3°  which  can  be  compared  directly 
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with  the  6.4°  maximum  rudder  angle  shown  in  Fiq.  20.  In  general,  the  optimal 
stochastic  controller  provides  a level  of  control  which  is  superior  to  that 
provided  by  a helmsman  given  by  Hooft's  model  as  implemented  in  eq.  (106)  and 
eq.  (107). 
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Figure  21.  Lateral  Offset  Response  to  Passing  Ship:  Helmsman  Control  at  H/T=' 


This  final  section  presents  the  results  of  an  evaluation  of  the  need  to 
adapt  path  control  systems  for  water  depth  dependent  changes  in  the  various 
terms  in  the  equations  of  motion.  In  general,  the  optimal  stochastic  path 
controllers  studied  here  provide  very  effective  control.  In  Section  2,  we 


showed  how  the  ten  coefficients  f . . and  y.  in  the  linear  equations  of 
motion  for  a ship  moving  along  a prescribed  straight-line  path  vary  with 


both  water  depth  and  vessel  speed.  In  Section  4,  we  studied  the  effect  of 
water  depth  and  vessel  speed  on  the  design  and  performance  of  the  optimal 
stochastic  path  controllers.  For  ships  such  as  the  Tokyo  Maru  which  are 
course  unstable  at  some  water  depth-to-draft  ratios  H/T,  we  concluded  that 
if  a design  is  to  be  optimal  at  a single  water  depth  (constant  C and  K 
gain  matrices)  it  should  be  designed  to  be  optimal  at  the  depth-to-draft  ratio 
at  which  the  ship  is  least  course  stable.  For  these  course -unstable  ships, 
however,  we  noted  the  need  for  an  adaptive  path  control  system  which  could 


automatically  adjust  the  C and  K gain  matrices  to  be  optimal  at  any 


operating  H/T.  In  these  systems,  the  coefficients  f^_. 


could  be 


identified  on-line  along  with  the  estimation  of  x . The  C and  K matrices 
could  then  be  updated  to  remain  optimal  and  reflect  the  changing  ship  charac- 
teristics with  changing  H/T. 

The  need  to  simultaneously  identify  all  ten  coefficients  in  eq.  (28);  i.e. 

^22'  ^23'  ^25'  ^21'  ^22'  ^32'  ^33'  ^35'  Y33'  1 32 ' cou-'-d  greatly  complicate  the 
design  of  an  adaptive  path  controller.  We  therefore  performed  a sensitivity 

study  to  see  how  important  it  would  be  to  adapt  for  changes  in  each  of  the 
individual  coefficients.  If  a specific  coefficient  were  not  to  be  identified 
on-line,  the  controller  design  would  most  likely  be  based  on  the  value  of 
that  coefficient  at  the  ship's  least  course-stable  water  depth-to-draft  ratio. 
The  study  was  therefore  conducted  by  first  using  the  OPTSYS  program  to  design 
optimal  controller  gains  C and  K for  the  Tokyo  Maru  using  the  Fe  matrix 
in  eq.  (58)  which  has  all  coefficients  f correct  for  H/T=°°  except  one. 

The  value  of  this  remaining  coefficient  was  set  at  its  value  -it  H/T=1.89, 
the  ship's  least  course-stable  condition.  This  simulates  the  gain  matrices 
which  a partially-adapted  path  controller  would  produce  if  it  were  to  identify 
all  coefficients  except  the  one  while  the  ship  is  operating  at  H/T=°°.  The 
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resulting  gain  matrices  were  then  used  again  in  the  OPTSYS  program  with  the 
Fe  matrix  which  was  fully  correct  for  H/T=°°  to  establish  the  resulting  RMS 
response  and  closed-loop  eigensystem.  This  is  the  response  of  the  ship  at 
H/T=°°  while  under  control  of  the  partially-adapted  controller  which  identifies 
all  coefficients  except  the  one.  This  process  was  repeated  for  each  of  the 
coefficients  in  succession  in  order  to  establish  a quantitative  sensitivity 
measurement  which  indicates  how  important  it  would  be  to  adapt  for  the  changes 


in  each  of  the  individual  coefficients.  An  adaptive  controller  which  would 
only  need  to  identify  and  adapt  for  changes  in  two  or  three  instead  of  ten 
coefficients  might  be  greatly  simplified  and  therefore  be  more  feasible. 


The  RMS  response  results  for  the  Tokyo  Maru  at  H/T=°°  with  the  partially- 
adapted  controllers  are  shown  in  Table  23.  Also  shown  for  comparison  are 
the  RMS  response  when  the  ship  is  controlled  by  the  optimal  controller  for 
the  existing  H/T=°°,  which  represents  the  fully-adapted  controller,  and  the 
RMS  response  with  the  optimal  controller  designed  for  the  ship's  least 
course-stable  water  depth  of  H/T=1.89,  which  represents  the  non-adapted 
controller.  The  results  are  arranged  in  the  order  of  increasing  RMS  cost  J; 
i.e.,  the  order  of  increasing  importance  for  adaptation.  Based  on  the  RMS 
cost,  seven  of  the  coefficients  (y31,  f35,  Y22,  f 33 , f22,  Y32-  and  f 32)  have 
very  little  effect  on  the  response  if  the  controller  does  not  adapt  for  their 
changes  individually.  Two  of  the  coefficients  y23  and  f^  have  a somewhat 
larger  effect  on  the  RMS  response.  The  final  coefficient  f , the  coefficient 
of  6'  in  the  r'  state  equation  (13),  has  by  far  the  greatest  importance 
based  on  this  measure.  Notice  that  the  three  most  important  coefficients  are 
all  in  the  r'  state  equation.  Notice  also  that  these  coefficients  do  not 
necessarily  experience  the  largest  numerical  variations  with  H/T  shown  in 
Table  4.  An  adaptive  control  system  could  probably  provide  most  of  the 
potential  improvement  in  performance  to  be  gained  with  adaptation  if  only 
f23<  ^25'  an<^  7^21  were  identified  on-line.  To  verify  this  conclusion,  we 
repeated  the  partially-adapted  controller  design  and  evaluation  process  for 
the  case  where  the  design  did  not  adapt  for  all  seven  of  the  coefficients 
Y(1#  f35<  Y22'  f33'  f22'  Y32'  and  f32  simultaneously.  This  represents  the 
situation  where  the  controller  would  adapt  only  for  changes  in  the  coefficients 


f 2 3 , f^,  anc*  Y2i  with  water  depth.  The  RMS  cost  J in  this  case  was  0.00683. 
The  cumulative  effect  of  not  adapting  for  the  changes  in  seven  of  the  coefficients 


fully-adapted 
optimal  at 
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all  adapted 
except 

all  adapted 
except 

all  adapted 
except  y22 

0.00965 

0.00966 

0.00967 

0.00960 

r 1 

0.00982 

0.00982 

0.00984 

0.00976 

6‘ 

0.00834 

0.00834 

0.00835 

0.00832 

n’ 

0.01607 

0.01609 

0.01607 

0.01610 

6 1 

0.04773 

0.04771 

0.04774 

0.04770 

<5  1 
c 

0.05406 

0.05404 

0.05405 

0.05405 

j 

0.00681 

0.00681 

0.00681 

0.00681 

all  adapted 

all  adapted 

all  adapted 

all  adapted 

except  fjj 

except 

except  Y32 

except  f 

32 

*• 

0.00946 

0.00966 

0.00955 

0.01032 

r ’ 

0.00968 

0.00963 

0.00968 

0.01047 

$• 

0.00831 

0.00831 

0.00830 

0.00850 

n’ 

0.01590 

0.01616 

0.01611 

0.01629 

6* 

0.04783 

0.04755 

0.04770 

0.04780 

6 1 

c 

0.05432 

0.05410 

0.05409 

0.05365 

J 

0.00681 

0.00682 

0.00682 

0.00683 

all  adapted 

all  adapted 

all  adapted 

non-adapted 

except  v 

except  f 

except  f 

optimal  at  H/T=1.89 

*• 

0.01029 

0.00834 

0.00846 

0.00850 

r’ 

0.01038 

0.00855 

0.00769 

0.00753 

0* 

0.00841 

0.00806 

0.00768 

0.00756 

n’ 

0.01775 

0.01539 

0.02066 

0.02238 

6* 

0.04690 

0.04824 

0.04556 

0.04533 

6 • 
c 

0.05224 

0.05599 

0.05279 

0.05249 

J 

0.00689 

0.00693 

0.00746 

0.00783 

Table  23.  RMS  Response  with  Partial ly-Adapte. . c’ontrollt  i at  H/T-  ■ 
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increases  the  cost  only  about  0.3  per  cent  and  is  no  worse  than  the  effect  of 
just  the  most  important  of  these  f ^ . 

The  closed-loop  eigenvalues  for  the  Tokyo  Maru  at  H/T=°»  with  the  partially- 
adapted  controllers  are  shown  in  Table  24.  Also  shown  for  comparison  are  the 
closed-loop  eigenvalues  when  the  ship  is  controlled  by  the  optimal  controller 
for  the  existing  H/T=°°  (the  fully-adapted  controller)  and  the  closed-loop 
eigenvalues  with  the  controller  designed  to  be  optimal  at  H/T=1.89  (the  non- 
adapted  controller).  The  coefficients  y21'  Y22'  Y31»  and  y32  do  not  affect 
the  feedback  control  gain  matrix  C and  thus  partially-adapted  controllers 
which  do  not  adapt  for  these  coefficients  have  the  same  closed-loop  eigenvalues 
as  the  fully-adapted  case.  These  coefficients  do  affect  the  filter  gains 
K and  therefore  affect  the  RMS  response  as  shown  in  Table  23.  The  partially- 
adapted  controller  closed-loop  eigenvalues  are  presented  in  Table  24  for  the 
remaining  coefficients.  These  results  are  presented  in  the  order  of  increasing 
importance  for  adaptation  based  on  RMS  cost  J as  shown  in  Table  23.  The 
dominant  closed-loop  eigenvalues  for  these  partially-adapted  controllers  and 
the  fully-adapted  and  non-adapted  controllers  are  shown  in  Fig.  23.  Again, 
coefficient  f is  the  most  important  single  coefficient  and  as  a minimum 
an  adaptive  system  should  identify  f on-line  and  update  the  control  and 
filter  gain  matrices  to  reflect  changes  in  f^3  with  water  depth.  Adaptation 
for  changes  in  coefficients  f and  y^  would  also  appear  justified. 
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Figure  23.  Dominant  Closed-loop  Eigenvalues  with  Partially-Adapted  Controllers 
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except  f 
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except  f 
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optimal  at 
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-0.8893 

-1.102 
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-6.997 
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-1.000 

-1.000 

-1.000 
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Table  24.  Closed-loop  Eigenvalues  with  Partially-Adapted  Controllers  at  H/T= 


7.  Conclusions 


The  principal  general  conclusions  of  this  work  are  summarized  here  for  reference. 

• With  sufficient  measurements,  the  optimal  stochastic  path  controllers 
as  developed  here  can  provide  very  effective  control  of  a surface 
ship  in  shallow  water.  Measurement  of  lateral  offset  is  necessary. 
Additional  measurements  of  heading  and  yaw  rate  are  needed  to 
provide  effective  control.  The  further  addition  of  a drift  angle 
measurement  from  a doppler  sonar  is  not  justified. 

• The  water  depth-to-draft  ratio  H/T  can  have  significant  effect  on 
the  performance  of  optimal  stochastic  path  controllers  when  they  are 
designed  for  a specific  water  depth  and  then  operated  over  the  feasible 
range  of  water  depths.  The  effect  of  speed  changes  on  the  performance 
of  fixed-gain  controllers  is  less  than  the  effect  of  changes  in  depth- 
to-draft  ratio. 

• The  effect  of  operating  an  optimal  stochastic  path  controller  at  a 
speed  which  is  slower  than  its  design  speed  appears  opposite  to  the 
effect  of  operating  in  a water  depth  which  is  shallower  than  its 
design  depth-to-draft  ratio.  As  a result,  operation  at  a condition 
which  is  slower  and  shallower  than  the  design  condition  provides 
performance  closer  to  optimal  than  operation  at  the  controller's 
design  speed  in  the  same,  shallower  depth  water. 

• Constant  gain  stochastic  path  controllers  should  be  designed  to  be 
optimal  at  the  ship's  least  course-stable  open- loop  water  depth-to- 
draft  ratio.  For  ships  which  are  course-stable  under  all  depth- 
to-draft  ratios,  this  approach  should  provide  a robust  design  which 
provides  very  effective  control  under  all  conditions. 

• For  ships  which  are  course-unstable  under  some  conditions,  constant 
gain  stochastic  path  controllers  which  are  designed  to  be  optimal  at 
the  ship's  least  stable  open-loop  condition  should  provide  effective 
control  under  all  conditions.  Significant  improvement  in  performance 
will  be  possible  for  these  ships,  however,  with  the  use  of  an  adaptive 
control  system  which  can  adjust  the  qains  to  remain  optimal  or  near 
optimal  for  all  operating  conditions. 
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• Two  features  of  the  control  and  filter  gains  could  be  utilized  to 
simplify  the  optimal  stochastic  controller  implementation  and  adapta- 
tion. The  control  gain  on  lateral  offset  is  constant  independent 
of  ship  type,  vessel  speed,  and  water  depth.  Filter  gains  for  the 
rudder  angle  estimate  are  zero  independent  of  ship  type,  vessel  speed, 
and  water  depth. 

• Digital  computer  simulation  confirms  the  very  effective  performance 
of  the  optimal  stochastic  path  controllers.  Simulation  also  confirms 
the  validity  of  using  first-order  shaping  filters  to  model  the 
design  passing  ship  yaw  moment  and  sway  force  disturbances.  The 
Kalman-Bucy  filter  is  very  effective  in  estimating  the  yaw  moment 
disturbance  but  much  less  effective  in  estimating  the  sway  force 
disturbance.  The  poor  sway  force  estimate  occurs  because  the  sway 
force  has  a relatively  small  influence  on  the  ship.  Fortunately, 
the  filter's  poor  estimate  of  the  sway  force  disturbance  does  not 
degrade  the  effectiveness  of  the  controller  for  exactly  the  same 
reason . 

• A comparison  of  simulations  of  the  Tokyo  Maru  controlled  by  an 
optimal  stochastic  path  controller  and  by  the  helmsman  heading 
control  model  presented  by  Hooft  and  adapted  here  for  path  control 
show  the  optimal  stochastic  path  controller  to  provide  a superior 
level  of  control.  This  conclusion  is  only  as  valid  as  the  human 
helmsman  modeling  used. 

• The  optimal  stochastic  path  controllers  developed  here  are  designed 
to  estimate  and  control  zero-mean  disturbances.  As  a result,  these 
designs  cannot  accommodate  a continuous,  non-zero-mean  disturbance 
such  as  a lateral  current.  These  disturbances  can  be  accommodated 
by  including  two  additional  shaping  filters  for  a constant  yaw  moment 
bias  and  a sway  force  bias  in  the  design.  We  are  proceeding  with  the 
development  and  study  of  these  controllers.  These  results  will  be 
reported  separately  later. 

” The  sensitivity  study  conducted  by  the  design  and  evaluation  of 
partially-adapted  controllers  shows  that  an  adaptive  path  control 
system  could  be  almost  as  effective  as  the  optimal  (99.7  percent  as 
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effective  for  the  Tokyo  Maru  based  on  the  RMS  cost  J)  if  it  were 
to  identify  and  adapt  for  changes  in  only  three  of  the  ten  coefficients 
in  the  system  equations  which  change  with  depth-to-draft  ratio. 


Coefficient  f of  in  the  r1  equation  is  of  greatest  impor- 


tance. Coefficient  f25  of  6*  and  y of  N'  in  the  r' 


equation  are  also  of  significance.  Changes  in  the  remaining  seven 
coefficients  with  depth-to-draft  ratio  have  a negligible  effect  on 
the  controller  performance.  We  are  proceeding  with  a development 
and  study  of  adaptive  ship  path  controllers.  These  results  will  be 
reported  separately  later. 
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